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PREFACE 


The  numerical  model  investigation  described  herein  was  authorized 
by  the  U.  S.  Army  Engineer  District,  Philadelphia,  on  19  February  1975. 
This  study  was  conducted  at  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  during  the  period  March  1975  to  January  1977  in  the  Wave 
Dynamics  Division  (WDD) , Hydraulics  Laboratory,  under  the  direction  of 
Mr.  H.  B.  Simmons,  Chief  of  the  Hydraulics  Laboratory,  and  Dr.  R.  W. 
Whalin,  Chief  of  the  Wave  Dynamics  Division. 

The  investigation  was  performed  and  this  report  prepared  by 
Mr.  H.  Lee  Butler,  WDD.  A significant  contribution  to  the  study  was 
made  by  Dr.  Donald  C.  Raney,  Professor  of  Engineering  Mechanics  at  the 
University  of  Alabama,  assigned  to  WES  under  terms  of  the  Intergovern- 
mental Personnel  Exchange  Act.  The  numerical  computations  associated 
with  this  work  were  performed  on  Control  Data  6600  and  7600  computers 
located  at  the  Los  Alamos  Scientific  Laboratory,  Los  Alamos,  New  Mexico, 
through  the  cooperation  of  C-Division. 

Directors  of  WES  during  the  course  of  the  investigation  and  the 
preparation  and  publication  of  this  report  were  COL  G.  H.  Hilt,  CE, 
and  COL  John  L.  Cannon,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 


CONTENTS 


Page 


PREFACE  1 

CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 

UNITS  OF  MEASUREMENT 4 

PART  I:  INTRODUCTION 5 

Background  5 

Objectives  8 

PART  II:  BASIC  EQUATIONS  OF  FLUID  FLOW 10 

Equations  of  Motion 10 

Numerical  Approach  12 

PART  III:  COMPUTATIONAL  TECHNIQUES  14 

Explicit  Finite  Difference  Solution  14 

Implicit  Finite  Difference  Solution  15 

Boundary  Conditions  16 

Numerical  Stability  18 

PART  IV:  COMPARISON  STUDY— MASONBORO  INLET  APPLICATION  21 

Method  of  Comparison  21 

Description  of  Computational  Area  and  Input  Data 21 

Discussion  of  Results  23 

PART  V:  CORSON-GREAT  EGG  COMPLEX 39 

Field  Surveys 39 

Computational  Grid 39 

PART  VI:  MODEL  VERIFICATION 43 

Corson  Inlet  43 

Great  Egg  Harbor  Inlet 53 

PART  VII:  ALTERNATE  IMPROVEMENT  PLAN  CONFIGURATIONS  62 

Inlet  Modifications 62 

Improvement  Plans  for  Corson  Inlet  63 

Improvement  Plans  for  Great  Egg  Harbor  Inlet  84 

PART  VIII:  CONCLUSIONS  Ill 

REFERENCES 113 

TABLES  1-5 

APPENDIX  A:  FORMULATION  OF  THE  WI  MODEL A1 


I 


2 


I ' 


CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (SI)  units  as  follows: 


Multiply 


feet 

miles  (U.  S.  statute) 

cubic  feet 

feet  per  second 

square  miles  (U.  S.  statute) 

square  feet  per  second 

acres 

feet  per  second  per  second 


By 


To  Obtain 


0.3048  metres 

1.609344  kilometres 
0.02831685  cubic  metres 
0.3048  metres  per  second 

2.589988  square  kilometres 

0.09290304  square  metres  per  second 

4046.856  square  metres 

0.3048  metres  per  second  per  second 


NUMERICAL  SIMULATION  OF  TIDAL  HYDRODYNAMICS 
GREAT  EGG  HARBOR  AND  CORSON  INLETS,  NEW  JERSEY 

PART  I:  INTRODUCTION 

Background 

1.  Great  Egg  Harbor  Inlet  and  Corson  Inlet  (GECI)  are  natural 
inlets  located  in  southeast  New  Jersey  (Figure  1)  between  Absecon  Island 
and  Ludlam  Beach.  The  inlets  are  located  in  areas  where  commercial  and 
recreational  fishing  as  well  as  resort  activities  are  an  important  part 
of  the  economy.  Both  inlet  systems  are  characterized  by  large  areas  of 
marshland  which  restrict  the  flow  to  a complex  system  of  channels.  The 
Corson  back  bay  area  contains  two  large  bodies  of  water,  Ludlam  Bay  and 
the  very  shallow  Corson  Sound.  The  southwest  portion  of  the  Great  Egg 
complex  is  characterized  primarily  by  two  large  water  bodies.  Great  Egg 
Harbor  Bay  and  Peck  Bay.  Great  Egg  Harbor  Inlet  system  joins  Absecon 
Inlet  system  in  the  north  and  Corson  Inlet  system  in  the  south  through 

a single  channel  at  Roosevelt  Bridge.  Some  freshwater  inflow  enters 
the  Great  Egg  system  from  the  west  through  Great  Egg  Harbor  Bay.  Corson 
Inlet  system  joins  Townsends  Inlet  system  in  the  south  through  a single 
channel  at  Sea  Isle  City. 

2.  Some  of  the  problems  which  affect  the  inlet  complex  are:  haz- 
ardous navigation  conditions,  inlet  instability,  shoal  conditions  within 
inlet  vicinity,  erosion  of  recreational  and  protective  beaches,  and 
tidal  flooding  of  oceanfront  property.  To  reduce  these  problems,  modi- 
fications to  Great  Egg  Harbor  and  Corson  Inlets  were  proposed  in  1966  by 
the  U.  S.  Army  Engineer  District,  Philadelphia  (NAPEN) , and  subsequently 
authorized  by  Congress  in  1970.  Detailed  development  of  plans  to  im- 
prove the  existing  conditions  at  GECI  was  deferred  until  completion  of 
this  investigation  to  study  the  tidal  hydrodynamics  of  existing  condi- 
tions as  well  as  plans  of  improvement  for  each  inlet.  Improvement  plans 
to  be  investigated  consist  of  the  installation  of  some  type  of  dual 
jetty  system,  including  a low-level  weir  on  one  jetty,  deposition  basin, 
and  channel  dredging. 
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, MARGATE  CITY 


Figure  1.  Location  of  GFCI  area  in  southeast  New  Jersey 


3.  This  study  was  funded  by  NAPEN  to  aid  in  evaluating  proposed 
modifications.  Plans  of  improvement  for  New  Jersey  inlets  and  beaches 
were  presented  by  NAPEN  in  a Navigation  and  Cooperative  Beach  Erosion 
Control  Study. ^ No  investigative  action  was  taken  until  a request  was 
made  of  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  in 
1973  to  provide  cost  and  time  estimates  for  a physical  model  study  of 
GECI.  A meeting  between  all  offices  concerned  was  held  to  fully  explore 
the  problem  and  issue  guidance  on  a feasible  study  plan.  As  a result 
NAPEN  requested  a numerical  model  investigation  of  the  GECI  area  in 
order  to  ascertain  tidal  hydrodynamic  evaluations  of  the  various  pro- 
posed plans. 

4.  The  first  major  problem  encountered  in  any  numerical  simulation 
study  is  the  selection  or  development  of  an  algorithm,  which  is  not  only 
capable  of  simulating  the  problem  under  study  but  can  also  be  economi- 
cally applied.  Most  numerical  simulations  of  tidal  hydrodynamics  asso- 
ciated with  inlets  and  bays  have  been  performed  using  explicit  finite 
difference  models.  An  explicit  scheme  is  one  in  which  the  solution  of 
the  difference  equation  at  a given  time  depends  solely  on  information 
computed  in  previous  time  steps. 

5.  The  first  computer  code  widely  recognized  for  treating  such 

problems  used  an  explicit  scheme.  The  code  was  developed  by  Reid  and 
2 

Bodine  (termed  the  RB  model)  and  was  used  extensively  in  simulations 
which  include  subgrid  topographic  features  and  possible  flooding  of  low- 
lying  areas.  An  extension  of  Reid  and  Bodine' s work  was  carried  out  by 

3 

Masch,  Brandes,  and  Reagan  to  help  evaluate  the  degree  to  which  mathe- 
matical models  can  be  used  to  predict  the  tidal  hydrodynamics  of  an  in- 
let system.  A suggestion  was  mad~  by  NAPEN  to  apply  Masch' s model  to 
the  GECI  area;  however,  past  experience  within  the  Wave  Dynamics  Divi- 
sion (WDD)  of  WES  had  shown  that  explicit  models,  in  many  instances, 
are  economically  infeasible  to  apply  to  real-world  engineering  problems 
covering  large  areas  due  to  stringent  stability  criteria  limiting  the 
calculation  time  step.  In  contrast  to  an  explicit  scheme,  an  implicit 
finite  difference  formulation  is  usually  not  restricted  by  an  associated 
stability  criterion.  For  an  implicit  formulation,  the  solution  at  a 


given  point  within  the  study  area  depends  upon  unknown  values  at  neigh- 
boring points  in  the  grid  system  as  well  as  upo*'  values  computed  in 
previous  time  steps. 

6.  Implicit  schemes  have  been  applied  by  WDD  to  a variety  of  prob- 
lems such  as  tsunami  propagation,  simplified  tidal  models,  landslide- 
generated water  waves,  and  storm  surge  calculations  summarized  by  Butler 

4 5 6 

and  Durham.  Implicit  formulations  were  used  by  Leendertse  ’ in  treat- 
ing problems  which  included  low-lying  areas  that  are  inundated  only  dur- 
ing a portion  of  the  tidal  cycle.  The  principle  reason  for  choosing  an 
implicit  formulation  is  one  of  economics.  A significantly  larger  time 
step  than  in  an  explicit  scheme  is  used  in  the  calculations  since  the 
time  and  spatial  increments  are  chosen  to  meet  required  accuracy  in  rep- 
resenting topographic  features  and  external  forcing  functions. 

7.  For  proper  application  of  a numerical  scheme  to  areas  as  large 
as  either  inlet  system  in  the  GECI  area,  an  implicit  model  was  deemed 
essential.  Since  there  were  no  implicit  codes  available  in  the  litera- 
ture which  could  simulate  flooding  of  low-lying  terrain  and  internal 
barriers,  such  as  dikes,  jettys,  weir  sections,  and  similar  features, 
the  development  of  such  a code  (WES  implicit  code  or  WI  model)  within 
the  scope  of  this  project  was  considered  necessary  and  was  included  as 
an  intrinsic  element  in  the  study.  In  order  to  evaluate  the  WI  model, 
a comparison  study  was  designed  to  check  its  reliability.  Both  the 

i Masch  version  of  the  RB  model  and  the  WI  model  were  applied  to  the 

Masonboro  Inlet  area  using  identical  grids,  bottom  topography,  and 
boundary  conditions.  Having  concluded  from  the  results  of  this  com- 
parison study  that  the  WI  model  was  reliable  and  economically  feasible 
to  apply  to  a large  inlet  system,  the  WI  code  was  selected  for  appli- 
cation to  the  GECI  area. 

Objectives 

8.  The  primary  objective  of  this  project  was  to  apply  a numerical 
model  to  Great  Egg  Harbor  and  Corson  Inlets,  N.  J.,  to  predict  quantita- 
tively the  hydrodynamics  (exclusive  of  salt  and  sediment  transport  and 


L 
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wave  action)  of  the  tidal  flow  through  each  inlet  for  various  improve- 


ment schemes.  Having  chosen  the  WI  model  for  application  to  GECI,  the 


following  tasks  were  accomplished: 


a..  The  WI  model  was  calibrated  for  existing  conditions  at 

each  inlet.  Tidal  elevations  and  current  velocities  were 
calculated  and  verified  by  comparison  with  prototype  data 
collected  by  the  staff  at  Rutgers  University. ^ 


Ij.  The  calibrated  model  was  applied  to  various  modified 

conditions  at  each  inlet  to  predict  tides  and  currents  at 
selected  stations.  Comparisons  are  made  to  indicate  the 
effects  of  each  improvement  plan. 


In  addition,  this  report  includes  a discussion  of  the  development  of  the 


WI  model  and  a presentation  of  the  comparison  study  made  for  the  Mason- 
boro  Inlet  application. 
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PART  II:  BASIC  EQUATIONS  OF  FLUID  FLOW 


Equations  of  Motion 

9.  The  hydrodynamic  equations  used  in  this  work  are  derived  from 
the  classical  three-dimensional  Navier-Stokes  equations  for  an  incompres- 

g 

sible  fluid.  The  details  of  the  derivation  can  be  found  in  Dronkers  or 
9 

Leendertse  and  are  not  repeated  here.  By  assuming  the  vertical  accel- 
erations are  small  and  the  fluid  is  well  mixed,  and  integrating  the  flow 
from  the  sea  bottom  to  the  water  surface,  the  normal  two-dimensional 
depth-averaged  form  of  the  equations  of  momentum  and  continuity  are  ob- 
tained. Expressed  in  a Cartesian  coordinate  system  (Figure  2)  the  in- 
tegrated equations  of  motion  are  given  as 


Figure  2.  Coordinate  system  for  problem  formulation 
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MOMENTUM 


3U  U 3U  V 3U  . . 3, 
3t  d 3x  d 3y  * 3x 


^(U2+V2)l/2  +F>< 


3V  U 3V  V 3V  ...  . 3 1) 

i7+7^  + 7^  + ,u+8di7 


^(u2+v2)1/2  + Fv 


CONTINUITY 


3,  3U  3V 

— - H + — • £ 

?t  3x  3y 


In  these  equations  U and  V are  the  vertically  integrated  transports 

per  unit  of  width  (flows)  at  time  t in  the  x and  y directions, 

respectively;  n is  the  water-surface  elevation  with  respect  to  the 

given  datum;  d = q - h is  the  water  depth  at  (x,y,t);  h is  the  land- 

surface  elevation  at  (x,y);  f is  the  Coriolis  parameter;  C is  the 

Chezy  frictional  coefficient;  F and  F are  terms  representing  ex- 

x y 

ternal  forcing  functions  such  as  wind  effects;  g is  the  acceleration 
due  to  gravity;  and  E,  is  a term  representing  the  rate  at  which  addi- 
tional water  is  introduced  into  or  taken  from  the  system  (for  example, 
through  rainfall  and  evaporation).*  The  approach  used  in  this  model 
differs  from  that  of  Leendertse  in  that  the  equations  of  motion  are 
formulated  in  terms  of  flows  rather  than  velocities.  In  addition,  the 
schemes  for  handling  the  flood-dry  procedures  and  subgrid  barriers  are 
unique.  It  should  also  be  pointed  out  that  the  convective  terms  of  the 

form  j -1^-  as  written  are  only  approximate.  A complete  derivation  of 

d 10 
the  equations  can  be  found  in  Wang  and  Connor.  The  complete  form  of 

the  convective  terms  is  significantly  more  difficult  to  treat  numeri- 
cally. Numerical  tests  were  carried  out  using  the  complete  form,  and 


* For  convenience,  symbols  and  unusual  abbreviations  ire  listed  and 
defined  in  the  Notation  (Appendix  D) . 
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negligible  differences  were  noted  for  the  types  of  applications  dis- 
cussed in  this  report. 

Numerical  Approach 

10.  The  basic  equations  of  motion  (1,  2,  and  3)  are  nonlinear 
partial  differential  equations.  A typical  study  region  possesses  a 
highly  complex  geometry  and  variable  boundary  conditions  exist  through- 
out the  system.  These  features  of  the  problem  make  a purely  analytical 
approach  intractable.  In  developing  a numerical  scheme  for  solving  the 
problem,  the  differential  expressions  are  replaced  with  finite  differ- 
ence expressions  to  operate  in  spatial  and  time  coordinates  on  specific 
points  of  a grid  system.  If  the  grid  size  approaches  zero,  the  solution 
of  the  difference  equations  must  converge  to  that  of  the  differential 
equations.  The  numerical  solution  schemes  presented  in  the  following 
section  satisfy  the  convergence  criteria.  Details  on  questions  of  con- 
vergence and  stability  can  be  found  in  references  cited  in  this  report. 

11.  Two  forms  of  finite  difference  solutions  to  the  equations  of 

fluid  motion  are  considered:  an  explicit  and  an  implicit  formulation. 
The  same  space-staggered  computational  grid  is  used  in  both  formula- 
tions. A rectilinear  mesh  is  placed  over  the  study  area  where  Ax  is 
the  height  of  a grid  cell  and  Ay  is  the  width.  Within  each  cell, 
the  following  assumptions  are  made:  (a)  the  value  of  q is  considered 
to  be  an  average  over  a grid  cell  centered  at  x = MAx  and  y = NAy 
and  at  time  t = kAt  ; (b)  the  value  of  M is  given  at  the  center  of 

the  lower  cell  face;  and  (c)  the  value  of  VN  M is  given  at  the  center 
of  the  right-hand  cell  face  (Figure  3).  In  addition,  the  water  depth 
and  frictional  coefficient  (Manning's  n or  Chezy  coefficient  C)  are 
also  defined  at  the  center  of  grid  cells. 

12.  The  following  notations  are  used  in  approximating  the  equa- 
tions of  motion  by  a system  of  difference  equations. 


Functional : 


UN  M = U(MAx,  NAy,  kAt);  Ax  = Ay  = A 
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<*>■ Am 


r’ 


D - FLOW/UNIT  WIDTH  IN 
X-DIRECTION  (U) 

A - FLOW/UNIT  WIDTH  IN 
Y- DIRECTION  (V) 

O - SURFACE  ELEVATION  („). 

WATER  DEPTH  (d), 

FRICTIONAL  COEFFICIENT  (C  or  n) 


Figure  3.  Cell  definition 


Time  Differential: 


9U  Uk+1/,2~  Uk"1//2 


at  (N,  M,  + 1/2) 


Spatial  Differential:  or  \ll)2 


Average  Terms: 


UN+1/2,M  “ UN,M+l/2  + UN-1,  M+l/2 


+ UN,M-l/2  + UN-l,M-l/2 


Angle  brackets,  / ^ or  ^ ^2  , are  used  to  indicate  that  terms 
maintained  in  differential  form  are  evaluated  with  centered  difference 
expressions  over  one  or  two  grid  cells,  respectively. 


PART  III:  COMPUTATIONAL  TECHNIQUES 


Explicit  Finite  Difference  Solution 

13.  The  explicit  solution  method  used  in  the  RB  model  is  a time- 

centered  difference  scheme  involving  a procedure  of  the  "leap  frog"  type 

for  computation  of  flow  and  water  levels.  Using  the  notations  given 

in  the  previous  section  and  applying  centered  differences  in  time  and 

1 

n can  be  expressed  explicitly  in  terms  of 


k+1/2  k+1/2 

space,  U , V 


functional  values  calculated  at  previous  time  levels.  An  abbreviated 
form  of  the  system  of  difference  equations  for  the  RB  code  can  be 
written  as: 


MOMENTUM 


= I 


Cfk-1/2  /-,|.\k-l/2 

,,‘♦1/1  . I M %-(^J  ♦ A«V‘-1/2 


- 4UC(^  * MFjj/c,,  AT  (N,  M * 1/2) 

=h'- 


iit»- 1/2 


'AfA'  AT  (N  t 1/2,  M) 


WHERE 


At  /BU\l,-,/2 


C 1 At/3U\ 

“ "‘Wt  T 2.a(d*//J 

W,  =j^(Uk‘l/2)*  + (Vk'1/2)*  I d* 

r , At/3V\k'1/2 

'C  + 4/3 

" \^n  2.21(d*)4/3  1 


I,  = (Uk,/2>1  ♦ (Vk'l/2)' 


1/2 


CONTINUITY 


*♦1  _ 
v * v 


-or- ■'] 


AT  (N,  M) 


14 


(4) 


(5) 


(6) 

(7) 

(8) 

(9) 

(10) 

(11) 


14.  The  calculations  for  each  time  step  are  divided  into  two 


halves:  the  flows  are  computed  during  the  first  half  of  the  time  step 
and  the  results  are  used  in  the  continuity  equation  to  calculate  the 
surface  elevations  during  the  second  half  of  the  time  step.  It  should 
be  emphasized  that  the  computed  flows  and  water  levels  are  determined 
one-half  time  step  apart,  yet  the  results  are  given  only  at  even  time 
levels. 

15.  The  equations  presented  above  are  applicable  for  cells  sur- 
rounded by  free  flowing  water  on  all  four  sides  (open  cells).  However, 
some  boundary  conditions  within  the  system  may  force  different  approxi- 
mations to  be  used  for  terms  such  as  the  advective  cross-product  terms 
9IJ 

(e.g.,  U~— ) . A detailed  discussion  of  the  explicit  scheme  can  be  found 
dx  2 

in  Masch,  et^  al . and  will  not  be  presented  here. 


16. 


Implicit  Finite  Difference  Solution 
To  solve  the  governing  equations  implicitly,  the  same  space- 


staggered  scheme  is  used.  The  implicit  code  employs  an  alternating- 
direction  technique  where  calculations  are  divided  into  two  parts.  The 
first  step,  or  1/2  cycle,  consists  in  solving  for  n and  U implic- 
itly; the  second  1/2  cycle  computes  q and  V implicitly.  The  omitted 
transport  in  each  1/2  cycle  is  assumed  constant  for  that  step.  Applying 
a centered  difference  operator  to  the  momentum  Equation  1,  and  the  con- 
tinuity Equation  3,  along  a grid  line  parallel  to  the  x-axis,  results 
in  a system  of  linear  algebraic  equations  whose  coefficient  matrix  is 
tridiagonal.  The  form  of  the  equations  for  the  first  1/2  cycle  is  given 
by: 
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(14) 


The  equations  for  the  second  1/2  cycle  are  similar  to  those  above  and 
are  not  presented. 

17.  In  this  numerical  scheme  the  surface  elevations  are  computed 
twice  within  each  complete  time  step;  the  flows  are  computed  once. 

Again,  these  equations  are  applicable  only  for  open  cells.  Various  ap- 
proximations must  be  introduced  where  boundaries  occur.  A complete  pre- 
sentation of  the  difference  equations  and  method  of  solution  is  given  in 
Appendix  A. 


Boundary  Conditions 

18.  Various  types  of  boundary  conditions  are  permissible  in  the 
present  system  of  computation  for  both  explicit  and  implicit  codes.  In 
both  codes,  boundary  conditions  at  the  tide  computation  boundary  (open 
boundary  conditions)  are  accomplished  by  setting  the  water  levels  nN  M 
as  prescribed  by  input  tables.  Flow  rates  may  be  specified  instead  of 
water  levels.  All  additional  boundary  conditions  relate  the  normal  com- 
ponent of  flow  at  the  boundary  to  the  state  of  the  water  level  at  the 
boundary.  These  conditions  may  be  grouped  into  two  classifications: 
water- land  boundaries  and  subgrid  barriers. 

Water- land  boundaries 

19.  Such  boundaries  are  prescribed  along  cell  faces,  hence  this 
condition  is  handled  by  specifying  U = 0 or  V = 0 for  those  cells 
when  impermeable  boundaries  exist.  In  estuarine  systems  with  large 
areas  of  low-lying  terrain  and  a significant  tidal  range,  many  areas 
alternately  dry  and  flood  within  each  tidal  cycle.  The  behavior  is 
simulated  by  making  the  location  of  the  land-water  boundary  a function 
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of  the  current  value  of  the  total  water  depth.  By  checking  the  water 
level  in  adjacent  cells  relative  to  the  ground  elevation,  a determina- 
tion is  made  as  to  the  possibility  of  inundation.  If  flooding  is  pos- 
sible, the  boundary  face  is  treated  as  open  and  computations  for  n , 

U , and  V are  made  for  that  cell.  The  drying  of  cells  is  simply  the 
inverse  process. 

Subgrid  barriers 

20.  Subgrid  barriers  are  defined  along  cell  faces  and  are  of  three 

types:  exposed,  submerged,  and  overtopping.  One  characteristic  of  such 

barriers  is  that  the  surface  elevation  is  computed  at  the  center  of  the 

cells  on  either  side  of  the  barriers.  A detailed  treatment  of  these 

3 

barriers  in  the  explicit  code  can  be  found  in  Masch,  et  al_.  The  fol- 
lowing discussion  is  limited  to  the  way  in  which  these  conditions  are 
simulated  in  the  implicit  code.  Exposed  barriers  are  handled  by  simply 
specifying  a no-flow  condition  across  the  cell  face.  This  type  of  bar- 
rier is  used  to  describe  dikes,  jetties,  and  similar  features  which  are 
impermeable  and  usually  of  a width  much  less  than  one  half  the  spatial 
grid  step.  Submerged  barriers  are  used  to  simulate  flows  across  such 
barriers  as  submerged  reefs,  dredged  material  disposal  banks,  pipelines, 
etc.  The  water  level  on  each  side  of  a submerged  barrier  must  always 
exceed  the  barrier  crest  elevation.  The  flow  over  a submerged  barrier 
can  be  controlled  in  a manner  similar  to  that  used  by  Masch  but  expe- 
rience has  shown  that  by  defining  a special  Chezy  coefficient  for  the 
barrier  face,  the  flow  over  the  barrier  can  be  simulated  without  intro- 
ducing unwanted  transients. 

21.  Overtopping  barrier  is  a terminology  used  to  distinguish  a 
barrier  which  can  be  submerged  during  one  portion  of  the  tidal  cycle  and 
totally  exposed  in  another.  Masch  used  a broadcrested  weir  formula  to 
describe  the  overtopping  nature  of  the  flow  and  then  the  submerged  weir 
formula  when  appropriate.  Since  a larger  time  step  is  used  in  the  im- 
plicit code  and  the  duration  time  of  overtopping  is  short,  a Chezy  for- 
mulation is  again  used  to  simulate  the  flow  across  the  barrier.  When 
the  barrier  is  exposed,  a very  small  Chezy  coefficient  (high  friction) 
is  used  to  "stop"  the  flow.  When  overtopping  occurs,  the  coefficient 
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is  increased  to  a specified  maximum  as  a function  of  the  water  level 
over  the  barrier.  As  the  water  level  decreases  at  a later  time,  the 
coefficient  is  decreased  accordingly. 

22.  The  approximations  necessary  to  implement  these  boundary 
conditions  for  the  implicit  code  are  presented  in  detail  in  Appendix  B. 

Numerical  Stability 

23.  For  an  explicit  solution  scheme,  the  grid  size  and  computa- 
tional time  step  are  related  through  a stability  criterion.  The  cri- 
terion associated  with  the  explicit  scheme  presented  here  is  given  by 
the  relation 


where  AS  is  the  mesh  size  and  d is  the  maximum  water  depth  occur- 

max  r 

ring  in  the  model.  This  approximate  condition  was  derived  from  expres- 
sions obtained  by  linearizing  the  problem.  When  the  nonlinear  terms  are 
included,  it  can  be  expected  that  the  time  step  will  require  further 
reduction. 

24.  Again  considering  the  linearized  equations,  it  can  be  shown 

that  the  implicit  difference  scheme  is  unconditionally  stable.  In  other 

words  the  space  and  time  steps  may  be  chosen  to  meet  required  accuracy 

in  representing  topographic  features  and  external  forcing  functions. 

8U  3U 

The  inclusion  of  the  nonlinear  advective  terms  (of  the  form  U ^ , V 
into  the  implicit  scheme  may  result  in  an  inherent  instability.  All  de- 
rivatives in  the  basic  equations  are  approximated  with  centered  differ- 
ences over  a single-grid  space  with  the  exception  of  the  advective  terms, 
which  are  computed  over  two  grid  cells.  Oscillations  of  the  water  level 
at  a grid  point  (period  4At)  may  occur  and  grow  unbounded. 

25.  A scheme  which  proved  effective  in  eliminating  these  instabil- 
ities was  the  use  of  a recursive  digital  filter  in  the  form 

k+1/2  'v  k+1/2  , k k-1/2 

n = an  + bp  + cn  (16) 
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where  n represents  the  latest  computed  water  level  and  n the  fil- 
tered water-level  value  actually  used  in  further  computations;  coeffi- 
cients a,  b,  c are  chosen  to  damp  oscillations  with  periods  on  the  order 
o.  4At  (corresponding  to  1/2  the  Nyquist  frequency),  while  permitting 
the  longer  period  wave  motion  to  remain  almost  undisturbed.  The  coef- 
ficients of  the  digital  filter  must  also  be  chosen  in  such  a way  as  to 
maintain  the  stability  of  the  filter.  By  comparing  the  results  of  ap- 
plying the  linearized  system  with  and  without  a filter  it  was  demon- 
strated that  filtering  has  a negligible  effect  on  the  solution  for  water 
levels  and  flows.  For  applications  presented  in  this  report,  values  of 
a = 0.6  , b = 0.3  , and  c = 0.1  were  selected  after  performing  a pa- 
rameter study  varying  the  coefficients  over  wide  range  from  (a,  b,  c)  = 
(0.89,  0.1,  0.01)  to  (0.4,  0.4,  0.2).  A more  detailed  discussion  of 
digital  filter  techniques  is  given  in  Appendix  C. 

26.  An  additional  instability,  which  is  termed  a "secondary  flow" 
phenomenon,  may  also  occur.  A discussion  of  this  problem  was  presented 
by  Vreugdenhil . 11  The  scheme  normally  employed  to  eliminate  this  in- 
stability is  the  inclusion  of  eddy-viscosity  terms  in  the  momentum  equa- 
tions of  the  form 


e 


(17) 


where  e is  taken  as  a constant  eddy-viscosity  coefficient.  Such  a 
term  is  actually  a representation  of  the  neglected  effective-stress 
terms  in  the  vertical  planes.  Effective  stresses  in  vertical  planes  not 
only  contain  viscous  and  turbulent  contributions  but  also  contributions 
from  the  momentum  flux  due  to  nonuniform  vertical  velocity  profiles. 

The  importance  of  these  terms  should  be  investigated  for  strongly  curved 
flows  with  considerable  divergence  or  convergence.  No  conclusive  theory 
has  been  formulated  as  to  the  mechanism  of  generation  of  secondary  flow 
and  its  relation  to  the  nonlinear  instabilities  involved.  One  method 
for  determining  whether  or  not  the  neglect  is  justified  is  that  of  hind- 
casting; that  is,  computations  are  made  with  and  without  these  terms  for 
a known  situation  (where  prototype  data  are  available).  The  value  of 
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e can  be  estimated  for  order  of  magnitude  from  the  relation 


e a 6Vs(U2  + V2)/C 


12 

Other  estimates  for  e have  been  given  by  Deardorff.  The  implementa- 
tion of  the  eddy-viscosity  terms  is  discussed  in  detail  in  Appendix  C. 


PART  IV:  COMPARISON  STUDY- -MASONBORO  INLET  APPLICATION 


I 


Method  of  Comparison 

27.  The  method  of  comparison  consists  of  applying  both  the  RB  and 
WI  codes  to  a 7.5-square-mile*  area  at  Masonboro  Inlet,  N.  C.  This  ap- 
plication was  selected  for  a number  of  reasons  including: 

a.  Good  prototype  data  were  available. 

1).  A physical  model  study  of  Masonboro  Inlet  was  performed 
at  WES  by  Sager  and  Seabergh^  making  available  additional 
data  for  calibration  and  verification. 

£.  Masonboro  Inlet  and  surrounding  area  is  a complex  system 
characterized  by  a combination  of  interconnected  channels 
together  with  a considerable  area  within  the  system  sub- 
ject to  flooding  above  mean  low  water. 

c[.  A jetty-weir  system  extends  seaward  from  the  northern 

outer  barrier  island  allowing  testing  of  the  subgrid  bar- 
rier features  available  in  both  codes. 

28.  The  conditions  (grid  area,  spatial  step,  boundary  conditions, 
etc.)  for  application  of  both  codes  were  identical.  Comparisons  of 
agreement  of  surface  elevations  and  depth-averaged  velocities  with  pro- 
totype measurements  were  made  for  each  scheme.  A comparison  was  also 
made  between  the  computer- generated  circulation  patterns  and  pictures 
of  surface  flow  patterns  taken  of  the  physical  model  (with  the  use  of 
confetti  particles).  This  application  differs  from  a coarse  grid-fine 
grid  approach  used  by  Masch  in  simulating  the  tidal  hydrodynamics  at 
Masonboro  Inlet. 

Description  of  Computational  Area  and  Input  Data 

29.  In  choosing  a computational  grid  for  simulating  the  hydro- 
dynamics at  Masonboro  Inlet,  a number  of  factors  were  considered.  Sea- 
ward open  boundaries,  at  which  the  ocean  tide  is  applied,  had  to  be  set 


* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (SI)  units  is  presented  on  page  4. 
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at  a distance  from  the  inlet  throat  area  such  that  transients  developing 
at  the  boundary  due  to  numerical  approximations  would  not  affect  the 
calculations  inland  from  the  entrance  channel.  Interior  forcing  bound- 
aries were  chosen  where  prototype  or  physical  model  data  exist.  The 
grid  size  had  to  be  chosen  small  enough  to  represent  the  topographic 
features  of  the  study  area  with  sufficient  resolution  to  permit  an  accu- 
rate simulation  of  the  tidal  flows.  A 300-ft  step  size  resulting  in  a 
grid  of  dimensions  41  x 57,  was  chosen  for  this  study  as  a compromise  be- 
tween a very  accurate  resolution  and  excessive  computational  costs.  A 
map  of  the  inlet  system  is  shown  in  Figure  4.  The  resulting  number  of 
computational  water  points  was  1721. 


Figure  4.  Areal  extent  of  com- 
putational grid  for  Masonboro 
Inlet,  N.  C. 


30.  The  jetty  system  protruding  from  the  outer  barrier  island  is 
composed  of  a weir  section,  with  a top  elevation  of  +2  ft  mlw  (Beaufort 
datum)  and  1000  ft  in  length,  extending  from  the  outer  island  to  near 
the  bend  in  the  jetty.  The  remainder  of  the  structure  is  impermeable. 
Exposed  barriers  were  used  to  represent  very  narrow  strips  of  high  land 
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in  the  system's  marsh  areas.  The  seaward  model  boundaries  were  set  at 
a distance  of  6,000  to  10,000  ft  from  the  inlet  throat. 

31.  Tidal  elevations  obtained  from  prototype  (survey  of  12  Sep  69) 
and  physical  model  data  were  imposed  at  boundary  lines  denoted  by  cir- 
cled numbers  in  Figure  4.  Tidal  elevations  at  boundary  3 within  the 
marsh  were  taken  as  approximately  the  same  as  those  in  the  innermost 
channels.  Bathymetric  surveys  of  the  computational  area  were  digitized. 
Frictional  coefficients  were  defined  by  assigning  number  codes  to  the 
various  types  of  terrain  and  applying  known  values  of  Manning's  n to 
these  number  codes.  The  program  uses  the  relation 


C-  1,49,1/6 


(19) 


to  calculate  the  Chezy  frictional  coefficient  as  a function  of  Manning's 
n and  total  water  depth.  Since  the  prototype  data  were  taken  on  a calm 
day,  no  wind  stress  was  applied  in  the  numerical  model.  Other  param- 
eters were  required  as  input  data,  including  such  items  as  gage  loca- 
tions, printout  and  plot  controls,  and  other  logic  control  parameters. 
Eddy-viscosity  terms  were  not  required  in  this  simulation  to  maintain 
numerical  stability.  The  selection  of  time-step  size  is  discussed  in 
the  next  section.  The  RB  and  WI  codes  use  essentially  the  same  data  but 
the  tabular  form  of  the  input  data  differs.  It  should  be  noted  that  a 
more  compact  form  for  the  input  data  was  programmed  for  the  WI  code. 


Discussion  of  Results 


32.  As  stated  previously,  both  the  WI  and  RB  codes  were  applied 
to  the  identical  computational  grid  approximating  Masonboro  Inlet  and 
associated  estuarine  system.  The  results  from  both  models  show  good 
agreement  with  prototype  measurements  at  all  gages.  The  locations  of 
velocity  ranges  and  tide  gages  for  the  12  Sep  69  prototype  survey  are 
given  in  Figure  5.  The  time-step  sizes  used  for  the  calculations  were 
3 sec  for  the  RB  model  and  90  sec  for  the  WI  model.  The  stability  cri- 
terion as  given  by  Equation  15  predicts  a 5.5-sec  time  step  for  stable 
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computations  involving  a linearized  set  of  equations.  Tests  were  made 
with  At  = 5 and  4 sec  and  instabilities  in  the  deeper  ocean  area  of  the 
grid  were  noted.  Calculations  with  At  = 3 sec  were  stable.  The  WI 
model  was  tested  for  At  = 45,  90,  and  180  sec;  the  results  for  the  two 
smaller  time  steps  were  essentially  the  same  and  the  results  for  At  = 
180  sec  showed  some  disparity.  Plots  of  the  current  velocity  at  a typ- 
ical gage  location  for  At  = 45,  90,  and  180  sec  are  given  in  Figure  6. 
On  the  basis  of  such  comparisons  a time  step  of  At  = 90  sec  was 
selected. 

33.  The  computations  were  begun  at  1350  EST  and  a fold-over  oc- 
curred at  2000  EST,  equating  input  boundary  conditions  at  this  point  in 
time  with  those  at  0730  EST.  This  procedure  may  cause  some  discrepancy 
in  the  results  but  is  required  since  the  models  must  "spin  up"  from  a 
quiescent  state  at  low  tide.  It  is  preferable  to  begin  the  computations 
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at  the  beginning  of  the  flood  phase  of  the  tidal  cycle.  The  four  tides 
used  to  drive  the  system  are  plotted  in  Figure  7.  Each  numbered  curve 
in  Figure  7 corresponds  to  the  appropriate  numbered  boundary  in  Figure  4. 
These  tides  were  obtained  from  either  prototype  records  or  physical 
model  data.  Results  of  the  surface  elevation  calculations  are  displayed 
in  Figure  8.  Both  models  describe  the  tide  gage  measurements  with  the 
same  degree  of  accuracy.  The  discrepancy  at  gage  2 in  the  ebb  phase  may 
have  resulted  from  two  problems.  First,  the  prototype  tide  gage  was 
located  within  a marina  behind  the  northern  barrier  island;  and  second, 
there  was  some  question  regarding  the  phasing  of  the  tidal  input  speci- 
fied at  boundary  2. 

34.  Figures  9-11  show  the  current  velocities  calculated  in  the 
numerical  models  relative  to  prototype  velocity  measurements  taken  at 
three  depths.  Prototype  velocity  data  were  taken  at  five  ranges  within 
the  inlet  system.  Velocities  have  been  plotted  as  being  either  flood 
or  ebb  according  to  the  general  direction  of  movement  across  each  range. 
Only  three  ranges  were  used  in  this  comparison  since  tidal  forcing 
boundaries  2 and  4 (Figure  4)  pass  through,  or  very  near  to,  velocity 
ranges  5 and  4 (Figure  5),  respectively.  Figure  12  is  included  to  de- 
pict results  calculated  at  stations  very  near  to  prototype  stations  on 
velocity  ranges  4 and  5.  It  also  is  noted  that  some  velocity  stations 
on  a range  are  separated  by  distances  less  than  a grid  step.  In  such 

• . cases  averages  of  flows  at  adjacent  cells  were  used  to  define  velocities 

at  the  gage  location.  Recalling  that  the  numerical  model  produces  a 
depth-averaged  unit  flow  rate,  integrated  velocities  were  obtained  by 
dividing  the  flow  by  the  local  water  depth.  Again,  good  agreement  was 
noted  for  both  numerical  models. 

35.  Figures  13-16  depict  a comparison  of  sample  circulation  pat- 
terns obtained  from  both  numerical  models  with  surface  current  patterns 
taken  from  the  physical  model  study  at  approximately  the  same  times. 

These  figures  were  constructed  by  plotting  the  instantaneous  depth- 
averaged  flows  at  each  grid  point  in  the  system  as  vectors.  The  di- 
rection and  length  of  each  vector  indicate  the  direction  and  magnitude 
of  the  flow.  Caution  must  be  taken  in  attempting  to  compare  the 
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Figure  9.  Comparison  of  current  velocity  agreement  with  prototype 
at  range  1 for  the  WI  and  RB  codes 
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Figure  10.  Comparison  of  current  velocity  agreement  with  prototype 
at  range  2 for  the  WI  and  RB  codes 
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Figure  11.  Comparison  of  current  velocity  agreement  with  prototype 
at  range  3 for  the  WI  and  RB  codes 
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Figure  13.  Comparison  of  flow  patterns  at  computational  hour  4 (1750  EST) 


Figure  14.  Comparison  of  flow  patterns  at  computational  hour  8 (0900  EST) 


patterns  at  computational  hour  12  (1300  EST) 


numerically  generated  flow  patterns  with  photographs  of  the  surface 
velocity  patterns  obtained  in  the  physical  model.  The  numerical  results 
portray  instantaneous  values  of  a vertically  averaged  flow,  whereas 
the  photographs  are  time  exposures  of  surface  currents  over  a period 
of  about  three  minutes  in  the  prototype.  For  Figures  15  and  16  the 
physical  model  photograph  times  differ  from  the  numerical  plot  times 
by  half  an  hour.  In  this  light  it  is  observed  that  the  major  features 
of  the  circulation  patterns  photographed  in  the  physical  model  are 
reproduced  in  the  numerical  models.  The  WI  model  appears  to  do  a 
better  job  in  describing  the  flow,  in  particular  the  flow  around  the 
jetty  and  over  the  weir  in  Figure  13,  channel  flows  in  Figure  14,  and 
the  direction  of  the  ebb  jet  in  Figure  15. 

36.  Table  1 presents  the  characteristics  and  run  statistics  for 
the  Masonboro  Inlet  simulation.  Recall  that  At  for  the  implicit 
scheme  is  the  time  for  a complete  computational  cycle,  that  is,  U and 
V are  computed  once  and  n twice.  In  the  explicit  scheme,  U , V , 
and  r|  are  computed  once  in  a time  step  of  3 sec.  The  relative  speed 
of  the  implicit  scheme  can  be  expressed  by  the  relation 
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(20) 


where  At^  and  Atg  are  the  time  steps  of  the  implicit  and  explicit 
schemes,  respectively,  and  is  the  execution  speed  of  the  implicit 

model  relative  to  the  explicit  model.  As  seen  in  Table  1,  a 15:1  ratio 
in  execution  time  was  obtained  for  the  time  steps  chosen. 

37.  In  the  application  of  numerical  models,  various  techniques 
have  been  used  in  previous  studies  to  "tune"  the  model  to  produce  a bet- 
ter agreement  with  prototype  data.  The  applications  presented  in  this 
report  have  been  made  in  a straightforward  manner.  Schemes  such  as 
placing  fictitious  submerged  barriers  within  channels  to  control  the 
flow  were  not  used,  nor  were  such  barriers  used  at  velocity  gage  lo- 
cations to  obtain  the  correct  flow  rate.  Higher  values  for  Manning's 
n were  used  to  represent  flow  restrictions  in  areas  where  the  actual 
land  formations  could  not  be  approximated  accurately  with  the  selected 
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grid  step.  Experience  with  this  study  indicated  that  by  correctly  ap- 
plying excitation  tidal  histories,  bathymetry,  and  standard  values  for 
frictional  coefficients,  both  of  the  numerical  models  (RB  and  WI)  could 
simulate  prototype  conditions  with  little  or  no  tuning.  If  tuning  is 
required,  it  can  be  accomplished  by  minor  adjustments  of  frictional  co- 
efficients. Results  of  the  comparison  study  clearly  demonstrate  the 


reliability  and  cost  effectiveness  of  tne  implicit  scheme  for  modeling 
tidal  hydrodynamics  of  a complicated  inlet-wetlands  system. 


PART  V:  CORSON-GREAT  EGG  COMPLEX 


Field  Surveys 

38.  A prototype  field  survey  at  GECI  was  performed  during  1974- 

7 

1975  by  the  staff  at  Rutgers  University,  New  Brunswick,  N.  J.  A com- 
plete description  of  the  geometric  characteristics  of  GECI,  prototype 
tide  data,  current  data,  and  hydrographic  survey  data  is  given  in  the 
Rutgers  report.  The  datum  for  all  data  was  taken  to  be  1929  SLD.  Tidal 
elevations  and  current  measurements  were  taken  at  each  inlet  over  a 
period  of  five  days,  interspersed  between  October  1974  and  June  1975. 
Ocean  tide  data  were  made  available  by  National  Ocean  Survey  (NOS) , 
an  agency  within  the  National  Oceanic  and  Atmospheric  Administration. 
Tides  in  the  study  area  are  semidiurnal  with  a mean  range  slightly 
greater  than  4 ft.  After  studying  the  available  data,  it  was  determined 
that  readings  taken  on  10  June  1975  at  Corson  Inlet  and  on  9 June  1975 
at  Great  Egg  Harbor  Inlet  were  of  higher  quality  than  those  taken  on 
other  days  and  should  be  used  for  model  verification.  In  addition,  the 
improvement  plans  were  modeled  using  the  same  forcing  tidal  functions  as 
those  used  for  verification,  providing  a base  for  comparison  relative  to 
existing  conditions.  A nonnegligible  wind  field  existed  during  the  June 
surveys  and  is  given  in  Table  2 for  both  Corson  and  Great  Egg  Inlets. 

As  applied  in  the  model,  the  wind  field  was  taken  as  time-dependent  but 
not  spatially  dependent. 

Computational  Grid 

39.  The  same  factors  considered  in  the  selection  of  the  Masonboro 
Inlet  grid  apply  to  the  GECI  application.  Due  to  the  size  of  the  GECI 
system,  prototype  data  were  not  collected  concurrently  for  both  inlets 
and  thus  it  was  necessary  to  model  the  inlets  separately.  As  a compro- 
mise between  resolution  accuracy  and  computational  effort,  a 300-ft  mesh 
size  was  selected.  This  resolution  allowed  a good  representation  of 
most  of  the  channels  within  each  system  and  was  also  adequate  for 
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describing  the  proposed  modifications.  Figures  17  and  18  show  the 
boundaries  of  the  model  limits  at  Corson  and  Great  Egg  Harbor  Inlets, 
respectively,  as  well  as  grid  alignment  and  pertinent  gage  locations. 
Most  of  the  islands  within  the  back  bay  regions  of  both  systems  are 
marshlands  and  hence  are  susceptible  to  flooding  during  high  tide.  The 
NOS  ocean  tide  data  were  recorded  at  the  Atlantic  City  Steel  Pier.  Grid 
dimensions  were  110  x 91  for  Corson  Inlet  and  152  x 106  for  Great  Egg 
Harbor  Inlet.  Those  portions  of  the  computational  grid  that  lie  seaward 
of  the  barier  islands  and  about  two  miles  north  and  south  of  the  inlets 
do  not  extend  to  the  limits  of  the  grid  dimensions,  thus  conserving 
the  computational  effort  (Figures  17  and  18).  The  bathymetry  was 
provided  by  the  Rutgers  survey. 

40.  Natural  boundaries  were  used  as  computational  boundaries. 
Typically,  these  boundaries  consist  of  nonflooding  barriers  through 
which  one  or  two  channels  flow.  In  the  field  survey,  tide  gages  were 
placed  in  the  channels  passing  through  these  boundaries  in  order  to 
obtain  proper  boundary  conditions.  For  Corson  Inlet,  the  Seashore  Line 
Railway  provided  a convenient  northern  boundary  and  the  county  road  from 
U.  S.  Highway  9 to  Sea  Isle  City  was  used  as  a southern  boundary.  The 
western  boundary  (nonflooding)  was  taken  near  the  Garden  State  Parkway. 
For  Great  Egg  Harbor  Inlet,  County  Highway  585  was  taken  as  a southern 
boundary,  Garden  State  Parkway  to  the  Somers  Point  tollgate  and  high 
land  from  there  to  Mill  Road  as  a western  boundary,  and  Mill  Road  to 
Margate  City  (County  Highway  563)  as  a northern  boundary.  Tide  gages 
were  located  at  Ludlam  Bridge  (Cor  2),  Seashore  Line  Bridge  (Cor  3), 
Roosevelt  Bridge  (GE  2),  Gateway  Marina  (GE  3),  Mainland  Marina  in  Dock 
Thoroughfare  (GE  4),  and  Alton  Marina  (GE  6).  Data  from  the  Mainland 
Marina  and  Alton  Marina  gages  were  averaged  to  produce  boundary  condi- 
tions at  Whirlpool  Channel  and  the  channel  between  Williams  and  Kiahs 
Islands  (GE  5).  As  previously  mentioned,  the  ocean  boundary  conditions 
were  taken  from  data  collected  from  the  ocean  tide  gage  at  Steel  Pier 
in  Atlantic  City,  N.  J.  (GE  1). 
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PART  VI:  MODEL  VERIFICATION 


Corson  Inlet 
" ■ - 


41.  The  objective  of  the  verification  phase  of  any  model  study  is 
to  demonstrate  the  model's  ability  to  produce  results  that  agree  with 
known  data  for  existing  conditions.  Such  a test  gives  confidence  for 
applying  the  model  as  a predictive  tool.  The  model  for  Corson  Inlet  was 
verified  by  simulating  existing  conditions  for  10  June  1975.  Compari- 
sons between  prototype  and  model  were  made  for  tidal  elevations  at  two 
gages  and  current  velocities  at  11  gages. 

42.  Numerous  problems  were  encountered  in  the  verification  pro- 
cess. The  source  of  most  problems  could  be  traced  to  the  prototype 
data.  A major  discrepancy  was  noted  in  phasing  of  the  tides  at  interior 
boundaries.  The  ocean  gage  tidal  data  supplied  by  NOS  agreed  closely 

with  predicted  amplitudes  and  phasing  given  in  the  National  Ocean  Survey 

14 

Tide  Tables,  East  Coast  of  North  America.  To  obtain  reasonable  re- 
sults it  was  necessary  to  shift  the  survey  tide  data  at  the  Sea  Isle 
City  and  Seashore  Line  Bridge  gages  in  time  such  that  they  agreed  with 
the  relative  phasing  given  in  the  NOS  tide  tables.  From  the  depth  data 
given  with  the  current  meter  readings  in  the  Rutgers  field  survey  re- 
port, it  was  apparent  that  the  survey  crew  had  great  difficulty  in  main- 
taining the  same  gage  location  throughout  the  survey  for  most  gages. 
Comparison  of  model  to  prototype  should  be  made  with  full  realization 
of  the  problems  enumerated  above. 

43.  Figure  19  is  a plot  of  the  three  tidal  hydrographs  used  as 
boundary  conditions  for  the  10  June  1975  verification.  An  initial  por- 
tion of  the  Seashore  Line  Bridge  data  was  missing  and  had  to  be  con- 
structed by  extrapolating  the  collected  data  from  0800  EST  back  to 
0400  EST.  Figures  20  and  21  show  the  comparison  of  tidal  elevations 
between  model  and  prototype  for  the  Strathmere  Coast  Guard  gage  and 
the  Whale  Creek  Marina  gage.  Figure  17  depicts  the  study  region  and 
locations  for  tidal  gages  and  velocity  stations.  The  discrepancy  in 
the  Strathmere  gage  results  was  apparently  due  to  an  error  in  relating 
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INPUT  TIDES  FOR  10  JUNE  I 97  5 - CORSON  INLET 
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Figure  19.  Forcing  tides  for  Corson  Inlet  model 


Figure  20.  Comparison  of  surface  elevation  agreement  with 
» prototype  for  the  Strathmere  tide  gage 


Figure  21.  Comparison  of  surface  elevation  agreement  with 
prototype  for  the  Whale  Creek  Marina  tide  gage 


the  prototype  tide  data  to  the  proper  benchmark.  This  gage  should  have 
an  amplitude  response  similar  to  the  ocean  gage  but  somewhat  lower  as 
indicated  in  the  tide  tables.  The  numerical  results  show  such  a be- 
havior for  the  Strathmere  gage. 

44.  Reproducing  current  velocities  is  a much  more  difficult  prob- 
lem than  reproducing  surface  elevations.  This  difficulty  arises  from 
several  sources:  (a)  the  numerical  model  calculates  a depth-averaged 
flow,  an  entity  which  cannot  be  measured  directly  in  the  prototype; 

(b)  prototype  velocities  cannot  be  measured  as  accurately  as  surface 
elevations  due  to  instrument  limitations;  and  (c)  the  numerical  model 
calculates  velocities  from  flow  data  computed  at  discrete  points  in  the 
grid  system  which  may  not  coincide  exactly  with  the  locations  of  proto- 
type velocity  gages.  Figures  22  and  23  depict  the  comparison  between 
prototype  and  calculated  velocities.  For  ease  of  inspection,  the  veloc- 
ities have  been  plotted  as  being  either  flood  or  ebb  according  to  the 
direction  of  movement  across  the  velocity  range.  Due  to  the  sparseness 
of  the  prototype  data  and  the  variation  in  peak  current  between  surface 
and  bottom  readings,  no  attempt  was  made  to  construct  smooth  curves 
through  the  data.  The  calculated  depth-averaged  velocities  are  repre- 
sented by  the  solid  dots  at  one-hour  intervals.  The  comparisons  show 
satisfactory  agreement  in  both  phase  and  magnitude.  Discrepancies,  such 
as  those  noted  for  the  inlet  throat  gage  in  Figure  22,  may  be  explained 
in  part  by  the  shifting  of  the  gage  location  in  the  field  survey. 

Gage  C in  Ben  Hands  Thoroughfare  was  the  only  gage  for  which  poor  agree- 
ment between  prototype  and  model  was  obtained.  The  poor  agreement  was 
apparently  due  to  a lack  of  sufficient  data  on  the  bathymetry  in  the 
back  bay  area,  in  particular,  the  data  relative  to  defining  the  proper 
flow  communication  between  Corson  Sound  and  the  Intracoastal  Waterway. 

45.  The  study  region  was  too  large  for  plotting  circulation  pat- 
terns on  a single  cathode  ray  tube  (CRT)  film  frame.  As  a compromise, 
an  area  consisting  of  3250  grid  points  around  the  inlet  throat  was 
plotted.  Figures  24-27  show  the  vector  flow  patterns  for  four  times 
in  the  tidal  cycle.  There  were  no  prototype  data  available  for  direct 
comparison  with  these  plots.  Two  types  of  circulation  patterns  are 


STRATHMERE  BRIDGE - GAGE  B 4T  OCEAN  DRIVE  BRIDGE,  NORTH  OF  INLET -GAGE  NC 


Figure  23.  Comparison  of  current  velocity  agreement  with  prototype  for  Corson  Inlet 

base  gages  B,  NC,  IN,  1W,  and  C 


Lgure  27.  Comparison  of  unit  flow  and  velocity  patterns  at  1900  EST  for  verification  runs  of  Corson  Inlet 


given:  (a)  instantaneous  depth-averaged  unit  flows  at  each  point  in  the 

plot  area,  and  (b)  instantaneous  depth-averaged  velocities.  The  flow 

patterns  given  in  part  (a)  of  these  figures  were  made  from  colored 

slides.  The  lightest  background  indicates  land  and  the  dark  background 

water.  The  shaded  area  indicates  terrain  which  was  exposed  at  low- 

water  level  but  was  flooded  at  the  time  appropriate  to  the  plot.  Such 

flow  plots  indicate  the  amount  and  direction  of  water  flowing  through 

a given  cross  section  at  a given  time.  A vector  with  length  equal  to 
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a cell  width  represents  a flow  of  50  ft  /sec.  The  velocity  patterns 
shown  in  each  figure  correspond  closely  to  surface  current  patterns, 
giving  a more  accurate  description  of  the  circulation  and  flushing 
characteristics  of  the  system.  These  plots  were  obtained  from  black/ 
white  microfilm  and  exposed  land  is  indicated  by  the  absence  of  a vec- 
tor. A vector  with  length  equal  to  a cell  width  represents  a velocity 
of  4 ft/sec.  Both  flow  and  velocity  plots  for  verification  at  every 
1/2  hour  of  the  computation  are  available  in  microfiche  form  and  are 
attached  to  the  back  cover  of  this  report  (Appendix  E,  Card  FI).  The 
examples  displayed  show  the  tendency  of  the  flow  to  remain  confined  to 
the  main  channels  in  both  flood  and  ebb  stage.  A major  characteristic 
of  the  flow  patterns  is  the  well-formed  eddy  behind  Peck  Beach  in  the 
vicinity  of  gage  NC. 

46.  Mean  tide  levels  and  volumetric  discharges  across  various  vel- 
ocity ranges  (Figure  17)  were  calculated  by  the  model  for  a 12-1/2-hour 
tidal  cycle  extending  from  0730  to  2000  hours  EST.  A tabulation  of 
these  results  for  10  ranges  is  given  in  Table  3.  Since  the  tidal  cycle 
is  asymmetric  (Figure  19)  the  discharges  during  flood  and  ebb  in  the 
above  table  do  not  agree.  Ranges  1,  2,  and  3 are  at  the  entrances  of 

the  back  bay  system.  The  net  discharge  of  water  into  the  back  bay  sys- 
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tern  was  approximately  3.0  x 10  ft  during  the  particular  cycle  under  in- 
vestigation. Ranges  4 and  5 indicate  that  the  total  inflow  through  the 
inlet  throat  during  the  cycle  under  investigation  is  split  between  the 
two  branches  as  follows:  70  percent  through  the  Strathmere  reach  and 
30  percent  through  the  Peck  Beach  reach.  No  prototype  data  were  avail- 
able to  verify  discharges  through  velocity  ranges. 
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Great  Egg  Harbor  Inlet 

47.  The  problems  encountered  with  the  Corson  Inlet  prototype  data 
were  prevalent  in  the  field  survey  data  for  Great  Egg  Harbor  Inlet 
(GEI) . The  model  for  GEI  was  verified  by  simulating  existing  conditions 
for  9 June  1975.  Comparisons  were  made  between  prototype  and  .nodel  for 
a single  tide  gage  and  10  velocity  gages.  The  study  region  for  GEI, 
showing  gage  locations,  is  depicted  in  Figure  18. 

48.  Again,  major  discrepancies  were  noted  in  the  phasing  of  the 
tides  used  at  the  interior  boundaries.  The  ocean  gage  tidal  data  sup- 
plied by  NOS  agreed  closely  with  the  predicted  amplitudes  and  phases. 

A phase  lead  existed  for  all  interior  prototype  gage  data  relative  to 
the  tide  tables.  Thus,  it  was  necessary  to  shift  the  tide  data  forward 
in  time  at  the  following  gages: 

Great  Egg  Phase  Lead 


Tide  Gage  No. 

Name 

min 

2 

Roosevelt  Bridge 

20 

3 

Gateway  Marina 

20 

4 

Mainland  Marina 

25 

6 

Alton  Marina 

10 

Figure  28  is  a plot  of  the  six  tidal  hydrographs  used  as  boundary  condi- 
tions for  the  9 June  1975  verification.  For  an  approximation  of  the 
tide  at  Whirlpool  Channel  (gage  5),  an  average  between  gages  4 and  6 was 
used  in  the  absence  of  prototype  data.  Figure  29  shows  the  comparison 
of  tidal  elevation  between  model  and  prototype  for  the  Ocean  City  Coast 
Guard  gage  (gage  7).  Again,  there  is  a discrepancy  in  the  gage  results 
which  is  apparently  due  to  an  error  in  the  benchmark  used.  The  ampli- 
tude response  indicated  in  the  tide  tables  agrees  with  the  numerical 
results.  Gage  7 provided  the  only  prototype  tide  data  for  use  in  veri- 
fying the  model. 

49.  Comparisons  between  prototype  and  calculated  velocities  are 
drawn  in  Figures  30  and  31.  The  comparisons  show  good  agreement  at 
most  gages.  More  recent  bathymetric  data  for  the  back  bay  area  would 
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Figure  30.  Comparison  of  current  velocity  agreement  with  prototype 
for  Great  Egg  Harbor  Inlet  base  gages  BW,  BE,  NC,  1W,  IN,  and  IS 


be  required  to  obtain  better  agreement.  This  is  particularly  true  for 
gages  S,  C,  N,  and  D in  the  Garden  State  Parkway  Bridge  areas.  The  ac- 
curacy of  the  field  survey  data  did  not  warrant  additional  effort  to 
obtain  a better  agreement. 

50.  Circulation  patterns  for  the  GEI  throat  area  (3760  points)  are 
presented  in  Figures  32-35  for  four  times  in  the  tidal  cycle.  Both  unit 
flow  and  velocity  vector  patterns  are  given.  The  flow  diagrams  were 
made  from  colored  slides  and  the  velocity  diagrams  from  microfilm.  Flow 
and  velocity  plots  for  verification  at  every  1/2  hour  of  the  computation 
are  available  in  microfiche  form  and  are  attached  to  the  back  cover  of 
this  report  (Appendix  E,  Card  F2) . The  examples  chosen  show  the  ten- 
dency of  the  flow  to  remain  confined  to  the  deep  entrance  channel  at 
flood  stage  and  a tendency  to  move  over  the  large  shoal  area  off  Ocean 
City  during  ebb  stage.  There  is  some  erratic  behavior  in  the  velocity 
field  (Figure  35)  in  the  channel  area  immediately  to  the  north  of  Ocean 
City.  Time  and  funds  were  not  sufficient  to  completely  resolve  the  rea- 
son for  this  behavior,  and  it  is  not  recommended  that  direct  quantita- 
tive comparisons  be  made  in  this  area  for  existing  conditions  and  the 
various  improvement  plans.  However,  it  was  ascertained  that  this  be- 
havior, whether  real  or  not,  had  a negligible  effect  on  computed  data 

in  the  inlet  throat  and  estuary.  No  prototype  data  were  available  for 
comparison  of  the  velocity  fields. 

51.  Mean  tide  levels  and  volumetric  discharges  across  various 
velocity  ranges  (Figure  18)  were  calculated  for  a 12-1/2-hour  tidal 
cycle  extending  from  077  EST  to  1930  EST.  A tabulation  of  results  for 
14  ranges  is  given  in  Table  4 (listed  as  base  condition).  As  for  Corson 
Inlet,  a check  was  made  on  the  net  discharge  of  water  into  the  system 
over  the  chosen  tidal  range.  Ranges  1-7  are  at  the  entrances  of  the 
back  bay  system.  The  net  discharge  of  water  into  the  system  was  ap- 

3 

proximately  2.7  x io  ft  during  the  particular  cycle  under  investiga- 
tion. Ranges  8-11  indicate  that  the  inflow  through  the  inlet  throat 
during  this  cycle  was  primarily  directed  toward  the  southern  portion 
of  the  system. 
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PART  VII:  ALTERNATE  IMPROVEMENT  PLAN  CONFIGURATIONS 


Inlet  Modifications 


52.  The  major  problem  conditions  affecting  the  GECI  area  are  pur- 
ported to  be: 

a.  The  presence  of  a shifting  sandbar  on  the  updrift  side  of 
Corson  Inlet,  which  has  almost  closed  the  inlet,  making 
passage  through  the  inlet  difficult,  hazardous,  and  at 
times  impossible. 

b_.  The  existing  channel  in  Great  Egg  Harbor  Inlet  is  unsafe 
for  navigation  except  under  ideal  conditions  of  tide  and 
wind. 


In  general,  the  outer  shoal  areas  adjacent  to  both  inlet  channels  are 
reported  to  produce  sharp  ground  swells  at  all  stages  of  tide.  Waves 
break  across  both  inlet  openings  during  low  stages  of  tide  which,  com- 
bined with  ground  swells,  make  passage  of  small  craft  almost  impossible. 
Many  local  boat  owners  have  relocated  elsewhere.  In  order  to  eliminate 
these  problems,  plans  of  improvement  for  each  inlet  have  been  consid- 
ered. Two  plans  for  Corson  Inlet  and  three  plans  for  Great  Egg  Harbor 
Inlet  were  recommended  for  testing  by  NAPEN. 

53.  The  basic  improvement  plans  authorized  for  each  inlet  are 
similar.  In  general,  they  provide  for: 

a.  A jetty  for  navigation  protection  and  beach  erosion  con- 
trol located  on  the  updrift  side  (northern  side)  of  each 
inlet.  The  jetty  includes  a low-water  weir  section  and  is 
designed  to  intercept  southerly  littoral  drift  and  prevent 
movement  into  the  inlet  channel  by  causing  the  drift  to  be 
trapped  in  a deposition  basin  for  later  transferral  across 
the  inlet  to  the  southern  beaches. 

b_.  The  construction  and  maintenance  of  specified  navigation 
channels. 

£.  A bulkhead  and  jetty  system  located  on  the  downdrift  side 
of  each  inlet. 

d.  An  optional  development  of  a landfill  area  adjacent  to  the 
bulkhead  for  public  recreational  use. 
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Improvement  Plans  for  Corson  Inlet 

54.  Tests  of  two  proposed  plans  of  improvement  were  conducted  with 
the  1974-1975  hydrographic  conditions  and  the  10  June  1975  tidal  condi- 
tions of  the  Corson  Inlet  system.  Figure  36  depicts  both  plans  which 
are  identified  as  Plan  AA  and  Plan  CC.  These  plans  included: 

a..  A 3225-ft-long  weir-jetty  system  extending  seaward  from 
Peck  Beach  with  a 2000-ft-long  low-weir  section  having 
a top  elevation  of  0.0  SLD. 

Ij.  A 33-acre  deposition  basin  (bottom  elevation  -29  ft  SLD) 
located  immediately  west  of  the  low-weir  section. 

c_.  A 300-ft-wide  entrance  channel  extending  to  the  -14  ft 
bottom  contour  and  maintained  at  -14  ft  SLD. 

d_.  A 100-ft-wide  bay  charnel  (to  be  maintained  at  -8  ft  SLD) 
linking  Strathmere  reach  and  the  entrance  channel. 

For  the  downdrift  side  of  the  inlet  the  improvement  plans  included: 

a.  For  Plan  AA,  a 3900-ft-long  bulkhead- jetty  system  extend- 
ing from  the  northern  end  of  Strathmere  along  the  inlet 
frontage  for  1200  ft  (bulkhead)  and  then  seaward  for 
2700  ft  (1100-ft  bulkhead  and  1600-ft  jetty) . 

t>.  For  Plan  AA,  the  fillet  between  the  jetty  and  the  bulk- 
head was  backfilled  with  sand  to  an  elevation  of  +6  ft 
SLD. 

c_.  For  Plan  CC,  the  bulkhead- jetty-fill  area  system  was  re- 
placed by  a 2800-ft-long  dogleg  jetty. 

Plan  AA  results 

55.  Figures  37  and  38  show  the  comparison  of  tidal  elevations  with 
Plan  AA  installed  using  the  verification  results  as  a base  condition. 
With  this  plan  installed  it  was  noted  that  tidal  elevations  were  not 
significantly  changed  with  the  exception  of  the  period  between  1100  and 
1600  hours  EST  (low  water).  During  this  period  the  flow  through  Middle 
Thoro  was  very  sensitive  to  any  change  in  input  data  such  as  a change 

in  the  tidal  boundary  conditions  or  the  introduction  of  an  improvement 
plan.  In  order  to  alleviate  the  problem,  more  accurate  hydrographic 
survey  data  for  the  back  bay  area  would  be  necessary.  In  particular. 
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Figure  36.  Proposed  improvement  plans  for  Corson  Inlet 
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Figure  37.  Comparison  of  surface  elevations  (plan  to  base)  for 
the  Strathmere  tide  gage 
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Figure  38.  Comparison  of  surface  elevations  (plan  to  base)  for 
Whale  Creek  Marina  tide  gage 


a better  description  of  the  flow  between  Corson  Sound  and  the  Intra- 
coastal Waterway  (which  includes  Middle  Thoro)  is  required. 

56.  Figures  39  and  40  depict  the  comparison  between  Plan  AA  and 
the  base  condition  for  current  velocities  at  base  gage  locations.  Gages 
located  well  inside  the  back  bay  system  displayed  no  significant  change 
in  velocity  outside  the  1100  to  1600  hours  EST  time  frame.  Gage  IN  lo- 
cated in  the  inlet  throat  showed  a 15  percent  increase  in  peak  flood 
velocity.  Gage  C located  in  Middle  Thoro  showed  a reversal  of  flow  be- 
tween 1200  and  1300  hours,  indicating  the  effect  of  the  numerical  prob- 
lem that  occurred.  Plan  AA  exhibited  a tendency  to  increase  the  veloc- 
ity of  the  flow  through  the  inlet.  Figure  41  depicts  the  comparison 
between  plan  and  base  conditions  for  current  velocities  at  a gage  lo- 
cated in  the  main  channel  of  the  jetty  system  (gage  E in  Figure  36). 

A second  gage  was  located  in  the  main  channel  midway  between  the  inlet 
throat  and  the  end  of  the  jetty  (gage  M in  Figure  36).  Although  the 
increased  flow  rates  cannot  be  quantified  due  to  the  numerical  problem 
occurring  between  1100  and  1600  hours,  the  tendency  to  increase  the  flow 
during  ebb  stage  is  apparent. 

57.  Figures  42-45  display  the  vector  unit  flow  patterns  for 
Plan  AA  at  the  same  times  in  the  tidal  cycle  used  for  displaying  the 
base  results  (Figures  24-27).  The  downcoast  jetty  in  the  velocity  pat- 
terns is  shown  extending  from  the  landfill  area.  The  prevalent  eddy 
behind  Peck  Beach  noted  in  the  base  results  was  still  maintained.  The 
plan  did  appear  to  keep  the  the  flow  well  confined  to  the  main  channel. 
The  flow  pattern  in  Figure  44  shows  the  development  of  a small  circula- 
tion cell  in  the  deposition  basin  at  the  end  of  ebb  stage.  Figures  42 
and  45  depict  a confusion  in  the  flow  pattern  at  flood  stage,  indicating 
the  mixing  of  a substantial  flow  coming  over  the  weir  with  the  main 
channel  flow.  A finer  mesh  size  in  the  throat  area  would  aid  in  allevi- 
ating the  computational  problems  in  this  area,  allowing  a more  accurate 
description  of  the  flow  pattern.  The  weir  action  appears  sufficient  to 
alleviate  sane  impoundment  on  the  updrift  side  of  the  weir- jetty  system. 
Both  flow  _nd  velocity  plots  for  Plan  AA  conditions  at  every  1/2  hour 

of  the  putation  are  available  in  microfiche  form  in  Appendix  E. 
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entrance  to  flat  creek -gage  A1  L entrance  to  flat  creek -gage  A 2 


Figure  39.  Comparison  of  current  velocities  (plan  to  base)  for 
Corson  Inlet  base  gages  Al,  A2,  Bl,  B2,  D1 , and  D2 
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Figure  40.  Comparison  of  current  velocities  (plan  to  base)  for 
Corson  Inlet  base  gages  B,  NC,  IN,  IW,  and  C 


58.  The  volumetric  discharges  across  three  of  the  ten  ranges 
listed  in  Table  3 are  plotted  in  Figure  46.  Again  the  numerical  problem 
is  evident  during  the  1100  to  1600  hours  EST  time  frame.  Outside  of 
this  time  interval,  the  net  discharge  is  essentially  the  same  with  and 
without  the  plan  installed,  indicating  that  the  plan  would  tend  to  have 
a negligible  effect  on  the  net  flow.  Due  to  the  numerical  difficulties 
experienced  in  running  the  plan,  the  time  integration  of  the  discharges 
would  give  suspicious  results  and  no  attempt  was  made  to  draw  informa- 
tion from  such  calculations. 

Plan  CC  results 

59.  Corson  Inlet  with  Plan  CC  installed  exhibited  behavior  similar 
to  the  Plan  AA  conditions.  Figures  37  and  38  show  the  comparison  of 
tidal  elevations  between  base  and  Plan  CC  conditions.  The  comments  on 
the  numerical  difficulties  given  in  previous  sections  apply  to  this  plan 
as  well.  Figures  39,  40,  and  41  give  the  current  velocity  comparisons 
for  Plan  CC.  Plan  CC  does  show  a tendency  to  produce  greater  current 
velocities  in  the  entrance  channel  during  ebb  stage  as  noted  in 

Figure  41. 

60.  Figures  47-50  display  the  vector  flow  patterns  for  Plan  CC  at 
four  times.  Without  the  landfill  of  Strathmere  Beach  (Plan  AA) , more 
flow  was  directed  during  flood  stage  over  the  shoal  area  between  the 
downcoast  jetty  leg  and  the  tidal  delta  located  just  inside  the  inlet 
(as  seen  in  Figures  47  and  50).  Both  flow  and  velocity  plots  for 

Plan  CC  conditions  at  every  1/2  hour  of  the  computation  are  available 
in  microfiche  form  in  Appendix  E.  Volumetric  discharges  are  plotted 
in  Figure  46.  Little  difference  is  noted  between  results  for  both 
improvement  plans. 

61.  Difference  plots  of  unit  flow  vectors  for  both  proposed  plans 
are  given  in  Figures  51-54.  The  vector  differences  were  formed  by  sub- 
tracting plan  results  from  base  calculations.  The  direction  of  the  re- 
sulting vector  will  be  reversed  relative  to  the  proper  flow  direction 
if  the  flow  magnitude  for  the  plan  is  greater  than  that  for  the  base 
condition.  The  effects  of  both  plans  are  noted  to  be  confined  to  a 
localized  area  near  the  improvement  structures.  During  peak  ebb  stage 
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Figure  46.  Comparison  of  volumetric  discharges 
(plan  to  base)  for  Corson  Inlet 
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Figure  54.  Comparison  of  unit  flow  difference  patterns  for  Corson  Inlet  Plans  M and  CC  at  1900  EST 


(Figure  53)  the  effects  appear  to  reach  into  the  Strathmere  channel  as 


well  as  into  Middle  Thoro.  These  effects  may  be  due  to  the  numerical 
difficulties  experienced  during  this  time  frame.  The  unit  flow  differ- 
ence plots  for  both  plans  at  1-hour  intervals  are  available  in  micro- 
fiche form  in  Appendix  E. 


Improvement  Plans  for  Great  Egg  Harbor  Inlet 


62.  Tests  of  three  proposed  inlet  modifications  were  conducted 
with  the  1974-1975  hydrographic  conditions  of  the  Great  Egg  Harbor 
Inlet  system  and  the  9 June  1975  tidal  conditions.  Figure  55  depicts 
Plans  A,  B,  and  Bl.  Plan  B1  consisted  of  only  the  upcoast  portion  of 

Plans  A and  B.  All  plans  tested  included  the  following  features:  j 

a_.  A 3000-ft-long  weir-jetty  system,  attached  to  an  existing 
jetty  at  Longport  and  extending  seaward  in  a southerly 
direction,  with  a 2200-ft-long  low-weir  section  having  a 
top  elevation  of  0.0  ft  SLD. 

b.  A 54-acre  deposition  basin  (bottom  elevation  -27  ft  SLD) 
located  immediately  west  of  the  low-weir  section. 

£.  A 300-ft-wide  entrance  channel  extending  seaward  to  the 
-14  ft  bottom  contour  and  bayward  to  Bascule  Bridge  at 
Ocean  City.  The  channel  is  to  be  maintained  at  -14  ft 
SLD. 

For  the  downdrift  side  of  the  inlet,  the  improvement  plans  included: 

a.  For  Plan  A,  a 5550-ft-long  bulkhead- jetty  system  extending 
from  the  northern  end  of  Ocean  City  along  the  inlet  front- 
age for  3450  ft  (bulkhead)  and  then  seaward  for  2100  ft 
(1100-ft  bulkhead  and  1000-ft  jetty). 

b. .  For  Plan  A,  the  fillet  between  the  jetty  and  bulkhead  was 

backfilled  with  sand  to  an  elevation  of  +6.0  ft  SLD. 

£.  For  Plan  B,  the  bulkhead-jetty-fill  area  system  was  re- 
placed by  a 4200-ft-long  dogleg  jetty. 

Plan  Bl  consisted  of  the  upcoast  structures  only  and  was  tested  solely 
for  investigatory  purposes. 

Plan  A results 

63.  Figure  56  shows  the  comparison  of  tidal  elevations  with  Plan  A 
installed,  using  the  verification  results  as  a base  condition.  With 
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Figure  56.  Comparison  of  surface  elevations  (plan  to  base)  for 
the  Ocean  City  tide  gage 


this  plan  installed  it  was  noted  that  tidal  elevations  were  not  signif- 
icantly changed.  The  plan  had  a tendency  to  reduce  the  tidal  elevation 
during  flood  stage  and  to  increase  it  during  ebb  stage  in  the  vicinity 
of  the  inlet  throat  but  the  tidal  amplitude  was  unchanged.  Mean  tide 
levels  for  various  locations  throughout  the  inlet  system  are  given  in 
Table  4,  and  it  can  be  seen  that  the  plan  caused  no  significant  mean 
tide  level  changes. 

64.  The  comparisons  between  Plan  A and  the  base  conditions  for 
current  velocities  at  base  gage  conditions  are  given  in  Figures  57  and 
58.  The  velocities  at  these  locations  were  changed  only  slightly. 

There  was  a general  tendency  for  maximum  ebb  and  maximum  flood  veloc- 
ities to  be  slightly  reduced.  The  base  gages  in  the  throat  area  are  not 
located  in  the  main  channel  and  thus  give  little  information  on  the  plan 
effects  in  this  area.  Figure  59  depicts  the  comparison  between  plan  and 
base  conditions  for  current  velocities  at  a gage  located  in  the  entrance 
channel  at  the  end  of  the  jetty  system  (gage  E in  Figure  55).  A second 
gage  was  located  in  main  channel  midway  between  the  inlet  throat  and  the 
end  of  the  jetty  (gage  M in  Figure  55) . Increases  in  maximum  velocity 
(30  to  70  percent)  are  noted  for  both  gages  during  flood  and  ebb  stages 
with  the  plan  installed. 

65.  Unit  flow  patterns  for  Plan  A at  the  same  time  in  the  tidal 
cycle  used  for  displaying  the  base  results  (Figures  32-35)  are  given  in 
Figures  60-63.  The  downcoast  jetty  in  the  velocity  patterns  is  shown 
extending  from  the  landfill  area.  Plan  A tended  to  keep  the  flow  con- 
fined to  the  main  channels  during  most  of  the  tidal  cycle.  During  ebb 
stage  the  flow  appeared  to  move  toward  the  deposition  basin.  Plots 
such  as  Figure  62  indicate  that  the  ebb  flow  might  well  erode  the  berm 
between  the  channel  and  the  deposition  basin.  The  flood  flow  over  the 
low-weir  section  during  higher  stages  of  the  tide  appears  substantial 
and  should  eliminate  any  impoundment  of  sand  on  the  updrift  side  of  the 
structure  by  allowing  the  sand  to  move  into  the  deposition  basin.  The 
downdrift  bulkhead  and  landfill  appeared  to  substantially  reduce  ebb 
flow  over  the  shoal  area  north  of  the  bulkhead.  Another  feature  of 
the  plan  was  the  formation  of  small  eddies  at  the  seaward  end  of  the 
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downcoast  jetty  during  ebb  stage  (Figure  61).  Both  flow  and  velocity 
plots  for  Plan  A conditions  at  ever  1/2  hour  of  the  computation  are 
available  in  microfiche  form  in  Appendix  E. 

66.  A tabulation  of  net  volumetric  discharge  across  various  veloc- 
ity ranges  is  given  in  Table  4.  As  for  the  similar  verification  results, 
the  calculations  summarized  in  Table  4 were  made  for  a 12-1/2-hour  cycle 
extending  from  0700  to  1930  hours  EST.  By  narrowing  the  inlet  entrance, 
the  net  flood  flow  through  the  inlet  was  reduced  by  9 percent.  The  ebb 
flow  was  likewise  reduced.  Time  histories  of  the  discharge  through 
range  1 (inlet  throat),  ranges  8 and  9 (Ocean  City  reach),  and  ranges  10 
and  11  (Longport  reach)  are  depicted  in  Figure  64.  The  plots  demonstrate 
the  reduction  in  volume  discharge. 

Plan  B results 

67.  Tidal  elevations  with  Plan  B installed  are  compared  with  base 
conditions  in  Figure  56.  Again  there  were  no  significant  differences 
between  plan  and  base  results  with  the  exception  of  a slight  phase  lag 
with  Plan  B installed.  Plan  B showed  less  change  than  Plan  A at  the 
Ocean  City  gage  but  somewhat  greater  differences  at  other  locations  as 
reflected  in  the  mean  tide  computations  tabulated  in  Table  4. 

68.  Current  velocity  comparisons  are  given  in  Figures  57-59. 

Plan  B gives  similar  results  to  those  noted  for  Plan  A for  the  verifi- 
cation gages.  This  can  be  expected  since  the  major  differences  should 
occur  in  the  landfill  area  which  was  removed  in  Plan  B.  Current  veloc- 
ities at  the  main  channel  gages  (Figure  59)  were  slightly  higher  during 
ebb  stage  for  Plan  B than  for  Plan  A. 

69.  Circulation  patterns  for  Plan  B at  four  selected  times  are 
displayed  in  Figures  65-68.  The  major  difference  between  this  plan  and 
Plan  A is  the  tendency  for  substantial  flow  to  seek  the  shorter  path 
across  the  large  shoal  area  north  of  the  downcoast  bulkhead.  Figure  65 
indicates  weaker  flow  streamline  than  corresponding  figures  for  Plan  A 
and  base  conditions  at  0900  hours  EST.  This  is  due  to  the  phase  lag 
noted  for  Plan  B.  The  marsh  area  north  of  Ocean  City  was  not  flooded 
at  1900  EST  because  of  the  phase  lag  (Figure  68).  During  the  first 
flood  stage  a numerical  problem  developed  in  this  shoal  region.  Since 
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the  simulation  for  each  plan  consisted  of  a multiple  number  of  runs  and 
successive  runs  used  restart  data  saved  in  the  previous  run,  a special 
damping  filter  was  applied  in  the  shoal  area  only  during  the  second 
flood  stage.  Thus,  the  flow  patterns  for  Plan  B during  hours  16  through 
18  EST  show  a smoother  pattern  in  the  shoal  area  than  that  noted  for 
hours  6 through  8 EST.  Both  flow  and  velocity  plots  for  Plan  B condi- 
tions at  every  1/2  hour  of  the  computation  are  available  in  microfiche 
form  in  Appendix  E.  The  ebb  flow  was  directed  toward  the  deposition 
basin  and  the  comments  on  Plan  A relative  to  this  problem  are  valid  for 
Plan  B as  well.  The  eddy  forming  at  the  seaward  end  of  the  downcoast 
jetty  was  more  prominent  for  Plan  B. 

70.  Table  4 also  includes  a tabulation  of  volumetric  discharge 
results  for  Plan  B.  The  net  flood  flow  through  the  inlet  was  slightly 
greater  than  that  for  Plan  A.  Plan  B had  a greater  effect  than  Plan  A 
in  altering  the  mean  tide  level  at  various  locations  through  the  system. 

The  mean  tide  was  slightly  reduced  throughout  the  backbay  system.  Time 
history  plots  of  discharge  are  given  in  Figure  64  for  those  ranges  used 
in  presenting  similar  Plan  A results.  These  results  demonstrate  that 

a small  phase  lag  was  experienced  for  Plan  B which  could  explain  why 
the  mean  tide  for  this  plan  was  lower. 

Plan  B1 

71.  Plan  B1  consisted  of  the  upcoast  structures  only  and  was 
tested  to  determine  if  the  flow  observed  in  other  plans  toward  the  de- 
position basin  during  ebb  stage  was  caused  by  the  downcoast  structures.. 

Plan  B1  per  se  is  not  a proposed  improvement  plan.  Tidal  elevations  at 
the  Ocean  City  gage  are  plotted  in  Figure  56  and  current  velocity  com- 


parisons at  the  verification  gages  are  given  in  Figures  57  and  58.  The 
results  are  very  similar  to  those  for  the  base  conditions.  What  this 
test  indicated  was  that  the  tidal  and  velocity  changes  observed  for 
Plans  A and  B were  caused  by  the  two- jetty  system  rather  than  by  the 
upcoast  jetty  alone.  The  velocity  increase  in  the  entrance  channel  area 
(Figure  59)  was  significantly  less  than  that  for  Plans  A and  B.  As 
shown  in  Table  4,  the  mean  tide  levels  for  Plan  B1  were  generally  closer 
to  base  conditions  than  were  those  for  Plans  A and  B.  The  same  was  true 
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for  discharges  across  the  various  ranges  throughout  the  model  (Fig- 
ure 64) . 

72.  Flow  patterns  for  Plan  B1  are  given  in  Figures  69-72  for  the 
same  times  in  the  tidal  cycle  used  for  displaying  results  from  other 
plans.  These  figures,  as  well  as  the  full  set  of  circulation  patterns 
available  in  microfiche  form  in  Appendix  E,  indicate  that  the  natural 
path  of  the  main  channel  funneled  the  ebb  flow  toward  the  upcoast  struc- 
ture. This  was  also  obvious  from  inspection  of  base  condition  flow  pat- 
terns (Figures  33  and  34)  relative  to  corresponding  flow  patterns  for 
Plans  A,  B,  and  Bl. 

73.  Difference  plots  of  unit  flow  vectors  for  the  three  plans 
tested  are  given  in  Figures  73-76.  As  for  Corson  Inlet,  the  differences 
were  determined  by  subtracting  plan  results  from  base  calculations.  In 
general,  the  effects  of  all  plans  are  noted  to  be  confined  to  the  local- 
ized area  near  the  improvement  structures.  The  effects  do  extend  rela- 
tively farther  into  the  back  bay  area  than  those  noted  for  the  Corson 
Inlet  modifications.  The  unit  flow  difference  plots  for  all  plans  at 
1-hour  intervals  are  available  in  microfiche  form  in  Appendix  E. 


Figure  74.  Comparison  of  unit  flow  difference  patterns  for 
Great  Egg  Inlet  Plans  A,  B,  and  B1  at  0900  EST 


Figure  75.  Comparison  of  unit  flow  difference  patterns  for 
Great  Egg  Inlet  Plans  A,  B,  and  B1  at  1100  EST 


Figure  76.  Comparison  of  unit  flow  difference  patterns  for 
Great  Egg  Inlet  Plans  A,  B,  and  B1  at  1900  EST 
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PART  VIII:  CONCLUSIONS 


74.  Both  Great  Egg  Harbor  and  Corson  Inlet  models  were  verified 
to  reproduce  satisfactorily  the  hydrodynamic  response  of  the  inlet  sys- 
tem to  specified  astronomical  tides.  The  verification  procedure  sub- 
stantiates the  model's  ability  to  reliably  predict  the  effect  of  the 
proposed  modifications  to  each  inlet  on  tidal  elevations,  current  ve- 
locities, and  circulation  within  the  study  area.  Although  numerical 
difficulties  occurred  in  the  back  bay  area  of  Corson  Inlet  during  a 
portion  of  the  tidal  cycle,  the  flow  patterns  and  hydrodynamic  data 
obtained  in  the  local  area  of  the  modifications  are  valid. 

75.  The  Corson  Inlet  model  results  indicate  that  relative  to  base 
conditions,  both  plans  caused  greater  ebb  flow  rates  in  the  entrance 
channel  while  maintaining  the  same  flood  flow  rates. 

76.  Greater  success  was  met  in  the  simulation  of  Great  Egg  Harbor 
Inlet.  The  concentration  of  ebb  flow  toward  the  upcoast  deposition  ba- 
sin and  structure  suggests  a potential  problem  with  the  proposed  plans. 

A separate  test  with  only  the  upcoast  elements  of  the  proposed  improve- 
ments installed  indicated  that  the  problem  was  associated  with  the  up- 
coast modifications  together  with  the  natural  alignment  of  the  main 
channel . 

77.  Now  that  the  models  are  formulated,  the  testing  of  additional 
plans  will  require  only  the  necessary  changes  in  the  code's  input  data. 
Future  testing  should  include  an  update  of  the  bathymetry  data  in  the 
back  bay  areas  of  GECI.  Such  testing  also  should  include  a parameter 
study  of  the  weir-deposition  basin  characteristics  once  a basic  plan 

is  selected.  The  results  presented  herein  are  sufficient  for  advance 
planning  purposes.  Detailed  design  studies  would  have  to  be  conducted 
separately,  including  a wave  refraction  analysis  to  assist  in  design 
for  final  jetty  alignment. 

78.  Results  of  the  comparison  application  at  Masonboro  Inlet, 

North  Carolina,  demonstrate  that  the  WI  model  can  reliably  and  econom- 
ically simulate  tidal  hydrodynamics  of  a complicated  inlet  system.  This 
conclusion  is  substantiated  by  successful  application  of  the  model  to 
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Great  Egg  Harbor  and  Corson  Inlets.  Statistics  on  computer  run  time  for 
the  application  of  the  WI  model  are  summarized  in  Table  5.  Success  of 
the  Corson  Inlet  application  was  limited  due  to  a lack  of  sufficient 
bathymetric  data  in  the  back  bay  area  to  resolve  the  proper  flow  com- 
munication between  some  of  the  interior  channels.  A few  numerical  prob- 
lems did  develop  in  some  of  the  plan  applications  and  various  filtering 
techniques  were  used  to  overcome  the  difficulties. 
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Dimensions 

Number  of  water  points 
Spatial  step  size,  ft 
Solution 

RB-Explicit 

Wl-Implicit 


41  x 57 
1721 
300 

At  Run  Time* 

sec  min 

3 60 

90  4 


* Based  on  18-hour  (prototype)  simulation  on  a 
CDC  7600. 


Table  2 

Wind  Field  During  Prototype  Survey 


EST 

Wind  Speed,  knots, 
and  Direction 

EST 

Wind  : 
and 

Speed,  knots. 
Direction 

Time 

6-9-75 

6-10-75 

Time 

6-9-75 

6-10-75 

0100 

E 5 

SE  6 

1300 

N 8 

SSE  11 

0200 

C* 

SSE  4 

1400 

NE  14 

SE  8 

0300 

C 

SSE  5 

1500 

EWE  15 

SE  8 

0400 

C 

SSE  4 

1600 

ENE  11 

S 10 

0500 

W 5 

C 

1700 

EWE  12 

S 8 

0600 

C 

C 

1800 

E 12 

SSW  9 

0700 

WSW  5 

E 6 

1900 

E 13 

SSW  10 

0800 

W 8 

E 6 

2000 

E 10 

S 10 

0900 

WNW  7 

ESE  6 

2100 

SE  9 

SSW  9 

1000 

NW  8 

E 8 

2200 

SE  11 

S 9 

1100 

WNW  10 

ESE  8 

2300 

E 8 

S 8 

1200 

NW  8 

SE  10 

2400 

SE  9 

SSW  9 

C = Calm 


Table  3 

Model  Mean  Tide  Levels  and  Discharges  at  Corson  Inlet  for 
10  June  1975  Between  0730  EST  and  2000  EST 


Range 

No. 

Mean  Tide 

Level,  ft,  SLD 

qin* ** 

ft3  x 106 

qout* 

ft3  x 106 

qtotal 

ft3  X 

1 

0.88 

872 

616 

256 

2 

0.84 

215 

166 

49 

3 

0.90 

20 

22 

-2 

4 

0.96 

553 

358 

195 

5 

0.95 

227 

182 

45 

6 

0.91 

112 

92 

20 

7 

0.91 

81 

65 

16 

8 

0.92 

193 

118 

75 

9 

0.90 

84 

57 

27 

10 

0.89 

123 

106 

17 

* For  ranges  1,  2,  and  3,  QJN  and  refer  to  flow  into  and  out 

of  the  back  bay  system;  for  all  other  ranges,  they  refer  to  flow 
under  flood  and  ebb  stage. 

**  A negative  sign  for  QT0'r^L  indicates  net  flow  out  °f  the  back  bay 


system  or  a net  ebb  flow. 


Model  Mean  Tide  Levels  and  Discharges  at  Great  Egg  Harbor  Inlet  (Base  and  Plan  Conditions) 


Table  5 

Computer  Run  Times  for  Tidal  Simulation  Code 


Grid 

Spatial 

Runtime/At 

Runtime/ 18 PH 

Inlet 

Dimensions 

Step,  ft 

sec* 

min** 

Masonboro 

41  x 57 

300 

0.167 

4.0 

Corson 

110  x 91 

300 

0.89 

21.4 

Great  Egg 

152  x 106 

300 

1.43 

34.2 

* All  computer  runs  were 
operating  system. 

executed  on  the 

LASL  CDC  7600 

using  the  CROS 

**  Run  time  per  18-hour  prototype  simulation. 


APPENDIX  A:  FORMULATION  OF  THE  WI  MODEL 


1.  To  form  the  finite  difference  approximations  of  the  differen- 
tial equations  given  in  Part  II  of  main  text,  a space-staggered  scheme 
is  used  wherein  flows  and  water  levels  are  defined  at  different  loca- 
tions in  the  grid  (Figure  Al) . The  technique  used  is  defined  as  a 
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Figure  Al.  Space-staggered  scheme 


multioperational,  alternating  direction  implicit  method.  The  first 
operation  (1/2  cycle)  computes  q and  U implicitly,  advancing  the 
time  from  kAt  to  (k  + l/2)At  . The  following  equations  for  the  first 
1/2  cycle  are  applied  columnwise  (along  lines  parallel  to  the  x-axis):* 


T - i!kJ*  \ 

Uk+!2  = yk-^  + At  f v - n>m'lM  n,m+3/2  n.m-^  / 

n,m+h  n,m+2i  L d*  2AS 

V /nk_!5  ..k-^  -v  ,*  r,  k+% 

2ASd*  ^n+l,m+*s  ^n-l.m+Jj  ® ^n.m+l 


* Additional  notation  used  within  this  appendix  is  defined  in  the 
text  and  is  not  included  in  the  list  of  general  notations  given  in 
Appendix  D. 
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..k-^ 

k+h  k-h  k-h  \ / o a c n,ra+4  + n.raA 
- n + n , - n )/2AS  + 2 • 

n,m  n,m+l  n,m  -z_2 

2(d*)^C 


((Uk"^  , )2  + V2)^}  + Fk 
n,m+V  x 


] 

k+h  _ k At  f.k+Js  ,.k+% 
^n.m  ^n.m  2As  jjn,m+%  " n,m-5; 
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where 


C = %(C  . + C) 

v n,m+l  n,m 


V = ^(Vk  j + Vk  , + Vk  , . + Vk  , ) 

n-^,m  n+h,m  n-^.m+l  n+^,m+l 

d*=  %(r)k  , + r|k  - h , - h ) 

n,m+l  n,m  n,m+l  n,m 


Gathering  terms  to  be  computed  along  line  n at  time  level  (k  + l/2)At, 
Equations  Al  and  A2  can  be  written  as 


k+h  — ,,k+h  k+h 

■at)  + a , U , + a , n , = B , 
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The  coefficients  are  defined  as  follows: 
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Assuming  that  nM  2 is  a given  water  level  at  the  lower  boundary  and 
k+^  M 

,2  is  a given  flow  at  the  upper  boundary,  the  matrix  form  of  the 


equations 

for 

line  n 

can  be  written 

as 

follows : 
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2.  The  system  of  equations  obtained  above  is  noted  to  be  tridiag- 
onal and  can  be  solved  with  a minimum  number  of  operations  by  defining 
the  process  of  elimination  as  a set  of  recursive  formulas.  Starting 
with  the  first  equation 


Uk+^  = -R  nk+5f  + S 
m+H  m m+1  m 


(A5) 


where 


r - • S = ■ 

m ' - ’ « - 

m+h  m+H 


and 
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Substitution  of  Equation  A5  into  the  second  equation  gives 
~am+%  ['Rm  nmfl  + Sm]  + %+! 
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Again,  the  flow  rate  can  be  expressed  as  a function  of  the  next  water 
level: 


,k+55  = -R 

m+3/2  m+1 


^1+2  + ^m+1 


(A7) 


where 


m+1 


am+2  p _ Bm+3/2+am+l^m+l 

* m+l=  - 

am+3/2+am+l  m+1  am+3/2+am+l  m+1 


In  general,  the  recursion  formulas  can  be  written  as: 


where 
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3.  In  the  FORTRAN  code,  fractional  subscripts  are  not  permitted. 
For  use  in  the  program  a new  index  system  is  adopted,  whereby  water 
levels,  flows,  etc.,  appearing  in  the  shaded  area  of  Figure  A1  have  the 
same  coordinate  index  (N,M).  Using  this  notation,  the  expanded  form  of 
the  recursion  formulas  for  the  first  1/2  cycle  can  be  expressed  as: 


A5 
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Coefficients 
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Solution 
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M M N,M 
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The  recursion  coefficients,  P , Q , R , and  S , can  be  computed  in 
succession  between  boundaries  on  line  N . Depending  upon  the  types  of 
boundary  conditions  that  are  applied,  various  approximations  of  these 
coefficients  are  required.  Having  calculated  the  set  of  coefficients, 
the  solution  equation  can  be  solved  for  all  surface  elevations  and 
flows  in  descending  order. 

4.  The  second  operation  (h  cycle)  computes  n and  V implicitly, 
advancing  the  time  from  (k+%)At  to  (k+l)At  . The  development  of  the 
recursion  formulas  for  this  \ cycle  is  very  similar  to  that  just  de- 
scribed. The  equations  for  the  second  cycle,  applied  row-wise  along 
a line  m , are: 
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5.  Using  the  same  notation  as  that  for  Equations  A10-A18,  the  re- 
cursion formulas  for  the  second  % cycle  can  be  written  in  expanded  form 
as: 

Coefficients 
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APPENDIX  B:  IMPLEMENTATION  OF  BOUNDARY 
CONDITIONS  FOR  THE  WI  MODEL 


1.  Boundary  conditions  used  in  the  WI  model  include  the  following 
types:  open,  closed,  barrier,  and  moving  boundaries.  To  implement 
these  conditions,  approximations  to  the  recursion  formulas  presented  in 
Appendix  A are  necessary.  Computations  are  made  along  segments  of  a 
grid  line  with  boundary  conditions  defined  at  both  ends  of  the  line  seg- 
ment. The  term  "lower  bound"  is  used  to  denote  that  end  of  the  computa- 
tional line  segment  which  corresponds  to  the  smaller  spatial  index  value. 
The  higher  indexed  line  segment  end  is  referred  to  as  the  upper  bound. 

Closed  Boundaries 


2.  In  the  approximation  of  an  impermeable  boundary,  it  is  assumed 
that  the  flow  normal  to  the  boundary  is  zero.  This  forces  the  assump- 
tion that  the  coastline  must  pass  along  cell  faces  where  flows  are 
defined.  At  a lower  boundary,  I , a closed  condition  is  computed  by 
setting  the  recursion  coefficients,  R and  S , to  zero,  namely. 


RI  = si  = 0 


At  an  upper  boundary,  I , closure  is  accomplished  by  setting 


Uj  = 0 or  Vj  = 0 


before  computing  the  set  of  solution  equations. 


Open  Boundaries 


3.  At  open  boundaries,  water  levels  are  given  as  functions  of  time 
(in  tabular  form),  thus  forcing  these  boundaries  to  pass  through  the 
center  of  the  cells  (where  the  water  level,  d = n - h , is  defined). 

To  implement  an  open  condition  at  a lower  bound,  for  a boundary  parallel 


to  the  y-axis  at  line  N and  row  M , the  following  approximations  are 
made  to  the  recursion  coefficients: 


P 


M 


0 


k+h  Fk 
nN,M  + Hj.M 


0 


Using  these  approximations,  and  SM  are  defined  accordingly.  Note 

that  the  computation  of  the  advective  term  would  require  values  of  U 
outside  the  computational  field.  This  problem  is  solved  by  using  a 
linear  approximation  of  the  differential  equation  in  such  a case.  The 
treatment  of  lower  open  boundaries  parallel  to  the  x-axis  is  similar. 

4.  At  an  upper  open  boundary  parallel  to  the  y-axis  at  (N,M),  the 
normal  flow  rate  at  the  first  interior  cell  is  approximated  and  written 
as 


U 


k+b 
N,M- 1 


= "m-1 


k+?5 

'n. 


+ s 


M-l 


(Bl) 


where  R and  S include  the  approximation 


0 


The  treatment  of  upper  open  boundaries  parallel  to  the  x-axis  is  similar. 


Barriers 


5.  Exposed  barriers  interior  to  the  system  are  impermeable  bound- 
aries and  are  treated  in  the  same  manner  as  that  indicated  for  closed 
boundaries. 

6.  Submerged  barriers  may  be  approximated  in  two  ways:  A Chezy 

formulation  or  weir  formulation.  The  weir  formulation  is  similar  to 

2 

that  used  by  Reid  and  Bodine.  The  flow  over  a submerged  barrier  can 


B2 


be  taken  as  that  for  a submerged  weir,  namely. 


QN-iW'l"!  (B2> 

where  is  the  flow  normal  to  the  barrier;*  is  the  depth  of 

water  over  the  crest  of  the  barrier;  and  H 2 are  the  water  levels 

on  the  two  sides  of  the  barrier  (both  of  which  exceed  the  barrier  crest 
elevation,  Z^);  and  Cg  is  an  appropriate  nondimensional  discharge  co- 
efficient. In  the  weir  formula,  is  taken  as 


“b  ‘ «H1  * H2>  - Zb 


(B3) 


and  the  sign  is  such  that  the  flow  is  directed  toward  the  low-head  side 
of  the  barrier.  In  the  computational  scheme,  the  flow  across  a sub- 
merged barrier  is  not  specified  directly;  instead,  a term  incorporating 
the  effect  of  weir  control  is  introduced  into  the  momentum  difference 
If  the  barrier  is  broad  and  parallel  to  the  y-axis  then 
Noting  that 


equation 

,2  ..2 
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gd 
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AS(CsDbr 
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from  inspection  of  B3,  the  frictional  effect  in  the  difference  equation 
is  replaced  by  the  right-hand  side  of  B4.  The  implementation  of  this 
scheme  is  accomplished  by  substituting 

.k-b. 


rd* 

AS 


u; 


N,M  1 


‘W 


for  the  frictional  term  in  D1  (see  Equation  A12). 

7.  The  Chezy  formulation  is  accomplished  by  defining  a special  co- 
efficient at  the  barrier  face.  When  the  appropriate  momentum  equation 
is  applied  at  the  designated  barrier  cell,  a special  value  for  the 

* Additional  notation  used  within  this  appendix  is  defined  in  the  text 
and  is  not  included  in  the  list  of  general  notations  given  in 
Appendix  D. 


B3 


I 


frictional  coefficient  (C  = Cw)  is  used.  In  this  scheme,  is  taken 

as  a function  of  the  water  depth  over  the  barrier  crest. 

8.  Overtopping  barriers  are  taken  as  those  which  can  dynamically 
change  from  a totally  exposed  barrier  to  a submerged  barrier  and  vice 
versa.  Due  to  the  larger  time  step  used  in  the  implicit  scheme,  the 
short  duration  of  overtopping  is  not  simulated.  Flags  are  set  within 
the  code  to  change  the  type  of  computation  at  the  barrier  face  when  the 
water  level  on  either  side  of  the  barrier  face  reaches  a designated 
height. 


Moving  Boundaries 

9.  Moving  boundaries  allowed  by  the  code  are  used  to  simulate  a 
flooding  coastline.  At  the  beginning  of  a simulation  run,  a small 
amount  of  water  (e^)  is  placed  on  areas  which  are  dry  initially,  but 
designated  as  possible  flood  cells.  When  the  water  level  at  a cell, 
adjacent  to  a "dry"  flood  cell,  is  rising  and  reaches  a height  greater 
than  (h  + e^) , the  closed  boundary  condition  defined  at  the  interface  is 
removed  and  the  appropriate  flow  quantity  is  computed.  For  a cell  to 
dry,  the  reverse  process  is  performed.  A cell  will  "dry  up"  if  the 
flows  through  all  four  cell  faces  are  zero.  It  is  assumed  that  when 
a cell  dries,  a thin  level  of  water  (of  height  e^)  is  allowed  to  re- 
main on  the  cell. 


APPENDIX  C:  METHODS  OF  ELIMINATING 
NONLINEAR  INSTABILITIES 


1.  Due  to  computational  considerations,  the  advective  terms  are 
not  approximated  completely  by  centered  differences  in  time.  Complete 
centered  differences  are  preferable  for  accuracy,  but  the  computational 
effort  would  be  excessive.  Thus,  a lower  order  approximation  is  taken, 
whereby  the  spatial  derivatives  in  the  advective  terms  are  evaluated  at 
an  earlier  time  level.  As  discussed  in  PART  III  of  this  report,  oscil- 
lations may  develop  due  to  these  lower  order  approximations.  Thus  far, 
there  has  been  no  conclusive  theory  regarding  the  physical  explanation 
of  these  instabilities.  In  the  development  of  the  WI  code,  two  schemes 
were  used  to  eliminate  this  problem. 

Digital  Filtering 


2.  As  mentioned  previously,  results  compiled  during  the  develop- 
ment of  the  code  indicated  that  the  advective  terms  introduce  a high- 
frequency  noise  of  period  4At  in  the  computation  of  the  surface  ele- 
vations. The  use  of  digital  filtering  techniques  was  suggested  by  the 
nature  of  the  problem,  namely,  the  construction  of  a "low-pass"  smooth- 
. ing  filter  to  damp  the  higher  frequencies.  For  ease  of  computation,  a 

second-order  recursive  filter  of  the  form 


+ bn^  + 


(Cl) 


was  adopted  as  discussed  in  PART  III.  A parameter  study  was  performed 
which  resulted  in  selection  of  values  for  coefficients  a , b , and 
c , * namely 

a = 0.6  , b = 0.3  , c = 0.1  (C2) 


* Additional  notation  used  within  this  appendix  is  defined  in  the 
text  and  is  not  included  in  the  list  of  general  notations  given  in 
Appendix  D. 


Cl 


The  associated  transfer  function  (reduction  factor  for  tidal  amplitude) 
for  this  selection  is 


H(f)  = 


a 


(C3 


where 


j 2irAtf 
z = eJ 


(C4 


The  filter  is  stable  according  to  the  criterion  given  in  Otnes  and 
Enochson;15  if  the  denominator  of  the  transfer  function  is  in  the  form 

z2  - hjZ  - h2  (C5 

and  the  point  h1  , h?  lies  in  the  triangle  shown  below  the  filter  is 
stable. 


It  also  can  be  shown  that  if  r is  a real  number  and 


k+p  p k 

n ^ = r*n 


(C 


with 


< 1 , 


(C 


the  filter  is  stable. 


where  f = 1/At  is  the  sampling  frequency.  Since  the  oscillations 
are  of  period  4At  , and  the  sampling  period  is  At  , then  the  reduction 
factor  at  f/f  = l/4At/l/At  = 0.25  is  0.526  percent.  The  effect  on  the 
amplitude  of  the  input  signal  (tide)  is  given  in  the  above  figure  for 
f/f  = 0.001  (corresponding  to  a 12. 5-hour  tidal  cycle)  to  be  a reduc- 
tion of  0.037  percent.  A similar  curve  for  the  phase  response  indicates 
an  effect  of  0.166  percent.  Shorter  period  components  of  the  signal  can 
be  introduced  by  folding  the  tidal  input  at  12.5  hours.  The  percentile 
effects  on  amplitude  and  phase  are  0.077  percent  and  0.35  percent,  re- 
spectively. Results  obtained  by  running  a linearized  system  with  and 
without  the  filter  show  no  discernible  differences  in  surface  elevations 
or  flows  throughout  a 12.5-hour  prototype  simulation. 

Eddy-Viscosity  Terms 

4.  The  inclusion  of  terms  to  abate  instabilities  caused  by  a 
"secondary  flow"  phenomenon  was  discussed  in  PART  III.  These  terms 


C3 


are  added  to  the  right-hand  side  of  the  momentum  equations  (see  PART  II, 
Equations  1 and  2)  and  may  be  written  as 
Equation  1 
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The  finite  difference  approximations  for  these  terms  are: 
Equation  1 
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Since  such  terms  involve  flows  only  at  earlier  time  levels,  they  appear 
in  the  computation  of  or  as  given  in  Appendix  A.  The  relation 


“ 6 V g(U2  +V2)/C 


(C 12) 


gives  an  approximation  of  the  value  of  e to  use  in  running  the  code. 
For  the  entrance  channel  area  at  Great  Egg  Inlet  where  instabilities 
occurred,  the  following  orders  of  magnitude  were  noted  for  the  param- 
eters appearing  in  Equation  Cl  2: 

■> 

g = 32.2  ft/sec“ 


C2  = 104  ft/sec2  -*■  c - S ft2/sec 


U,V  = 100  ft  sec 


C4 


Computer  runs  were  made  varying  e (over  a range  of  2 to  20) , and  the 
results  were  noticeably  affected  in  both  amplitude  and  phase  as  e was 
increased  beyond  e = 5 , It  is  important  to  choose  the  smallest  value 
of  e such  that  the  instability  problem  is  cured  with  minimal  effect 
on  the  accuracy  of  the  computations. 


APPENDIX  D:  NOTATION 
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d 


max 

f 

F , F 
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k 


M,N 
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n ,n 
x y 

t 

u,v 

x,y,z 
Ax,  Ay 


AS 

At 
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n 
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Chezy  frictional  coefficient 
Total  water  depth 

Maximum  water  depth  at  any  location  and  time  in  the  system 

Coriolis  parameter 

External  forcing  functions 

Acceleration  due  to  gravity 

Land-surface  elevation 

Time  increment  counter 

Indices  denoting  spatial  increments  in  the  x and  y direction 
Manning's  frictional  coefficient 
Averaged  Manning's  coefficient 
Time 

Vertically  integrated  horizontal  transport  (flow)  components 
Cartesian  coordinates 
Spatial  increment 

General  spatial  increment  (grid  size) 

Time  increment 
Eddy-viscosity  coefficient 

Water-surface  elevation  with  respect  to  given  datum 
Rainfall  minus  evaporation 
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